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Abstract 



We consider quite general /i-pseudodifferential operators on R™ with small 
random perturbations and show that in the limit h — * the eigenvalues are 
distributed according to a Weyl law with a probabality that tends to 1. The 
first author has previously obtained a similar result in dimension 1. Our class 
of perturbations is different. 



Nous considerons des operateurs /i-pseudodifferentiels assez generaux et 
nous montrons que dans la limite h — > 0, les valeurs propres se distribuent 
selon une loi de Weyl, avec une probability qui tend vers 1. Le premier 
auteur a deja obtenu un resultat semblable en dimension 1. Notre classe 
de perturbations est differente. 
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Introduction 





This work can be viewed as a continuation of [7], where one of us studied random 
perturbations of non-self adjoint /i-pseudodifferential operators on R and showed that 
Weyl asymptotics holds with a probability that is very close to 1. In the present 
work we consider the multidimensional case and weaken some of the assumptions 
in [7] (like independence of the differentials and analyticity of the symbol). Our 
random perturbations are slighly different however, in |7j they are given by a random 
potential while here they are rather given by a random integral operator. 

Before continuing the general discussion, we fix the framework more in detail. 
We will work in the semi-classical limit on R n . Write p = (x, £) and let m > 1 be 
an order function on the phase space R 2 '^: 

3C > 1, iVo > such that m{p) < C (p - p) N °m(p), (1.1) 
V P , fi e R 2n , ( P -p) = ^l + \ P -p\ 2 . 

The corresponding symbol space (cf [2]) is then 

S(R 2n ,m) = {ae C°°(R 2n ); \d a p a(p)\ < C a m(p), p e R 2 ™, a e N 2 ™}. (1.2) 
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Let 

P{p; h) ~ p(p) + h Pl {p) + ... in S(R 2n , m). (1.3) 
Assume 3 z G C, C > such that 

\p(p)-z \>m(p)/C , PGR 2 ". (1.4) 

Let £ denote the closure of p(R 2n ) so that £ = p(R 2n ) U £oo, where Soo C C is the 
set of accumulation points of p in the limit (a;,£) = oo. 

For h > small enough, we also let P denote the ^i-Weyl quantization, 

Pu(x) = P w (x,hD x ,h)u(x) = ^— J j e^ x ~yyp(^-, V ;h)u(y)dyd V . 

Let Q d C \ Sqo be open simply connected and not entirely contained in E. 
Then, as we shall see, 

1° <j(P) fl Q is discrete for h > small enough, 
2° Ve > 0, 3 h(e) > 0, such that 

a{P) nQcEi D(0, e), < h < h(e). 

Here D(0, e) denotes the open disc in C with center and radius e > and we equip 
the operator P with the domain H(m) := (P — Zoj^^L^R")), where the operator 
to the right is the pseudo differential inverse of P — z (see [2] and [7j). 
If P is selfadjoint (so that p is real-valued) we have Weyl asymptotics: 
For every interval I C £1 with vo1r,2, 1 (p~ 1 (9/)) = 0, the number N(P,I) of 
eigenvalues of P in / satisfies 

N(P, I) = T^w (vol ip-\l)) + o(l)), h - 0. (1.5) 

This result has been proved with increasing generality and precision by J. Chazarain, 
B. Helffer-D. Robert, and V. Ivrii. (We here follow the presentation of j2] where 
references to original works can be found. The corresponding developement for 
selfadjoint partial differential operators in the high energy limit has a long and rich 
history starting with the work of H. Weyl (14j.) A very simple and explicit example 
is given by the harmonic oscillator P = \{{hD) 2 + x 2 ) : L 2 (R) L 2 (R), P(x,£) = 
p(x,^) = |(x 2 + £ 2 ). In this case the eigenvalues are given by \k{h) — (k + \)h^ 
k = 0, 1,2, ... 

In the non-self adjoint case, Weyl asymptotics does not always hold. If P is 
a differential operator with analytic coefficients on the real line, then often the 
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spectrum is determined by action integrals over complex cycles, having nothing to 
do with volumes of subsets of real phase space. A simple example of this is given 
by the non-self adjoint harmonic oscillator, 

P = ^((hD) 2 + ix 2 ) : L 2 (R) -> L 2 (R), (1.6) 

whose spectrum is equal to {e i7T ^(k + ^)h; k G N}; This is easy to see by the method 
of complex scaling, or by applying the general multidimensional result of [11J. In 
this case, we have Eoo = 0, and E is the closed 1st quadrant. Clearly the number of 
eigenvalues in an open set T <<= C intersecting the 1st quadrant, whose closure avoids 
the ray given axgz = j is equal to zero while the corresponding Weyl coefficient 
vol(p _1 (r)) is not. (Further results about the non- self adjoint harmonic oscillator 
have been obtained by E.B. Davies and L. Boulton, see [T| and further references 
given there). 

However, in this case and for quite a general class of /i-pseudodifferential opera- 
tors in one dimension, it was shown by one of us in jTj that if we replace the operator 
P by P + 5Q U , where < 8 <C 1 varies in a suitable parameter range and is a 
random potential of a suitable type then we do have Weyl asymptotic in the interior 
of E with a probability that is close to 1. The book j3] of M. Embree and L.N. 
Trefethen as well as the paper [13J by L.N. Trefethen and S.J. Chapman contain (in 
our opinion) numerical examples where one can see the onset of Weyl-asymptotics 
after adding small random perturbations. 

In this work we establish similar results in arbitrary dimension that we shall now 
describe. Let < fh, m < 1 be square integrable order functions on R 2n such that 
m or m is integrable, and let S G S(m), S G S(m) be elliptic symbols. We use 
the same symbols to denote the /i-Weyl quantizations. The operators S, S are then 
Hilbert-Schmidt with 

||£||hs,||£||hs-/^, 

where ^ indicates same order of magnitude. Let ei, e^, and ei, ejj, ... be orthonor- 
mal bases for L 2 (R n ). Our random perturbation will be 

where a^k are independent complex 7V(0, 1) random variables, and e~je* k u = (u\ek)ej, 
u G L 2 . In the appendix Section [13] we show that up to a change of the set of 
independent jV(0, l)-laws, the representation (18. II) is independent of the choice of 
bases e} and e.j. 
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Let 

M = C 1 h' n , (1.8) 

for some C\ ^> 1. Then, as we shall see in Section [HI we have the following estimate 
on the probability that Q be large in Hilbert-Schmidt norm: 

P{\\Q\\ls > M 2 ) < Cexp (-h- 2n /C), (1.9) 

for some new constant C > 0. In the following discussion we may restrict the atten- 
tion to the case when ||<3u||hs < M. We wish to study the eigenvalue distribution 
of P + 5Qu for 5 in a suitable range. 

Let r (<= Q be open with C 2 boundary and assume that for every z G dT: 

S 2 := p^ 1 {z) is a smooth sub-manifold of T*R™ on (1-10) 
which dp, dp are linearly independent at every point. 

The following result will be established in Section [TUJ 

Theorem 1.1 Let V d Q be open with C 2 boundary and make the assumption 
U.10\) . For 0</i<l, let 5 > be a small parameter such that 

< 5 < h 3n+1/2 . 

For some small parameter < e < 1 assume < f < 1 (or equivalently 

5 > e~ e ^ DhS > for some D ^> 1, implying also that e > Chlnh 1 for some C > 0). 
Then there is a constant C > (that is independent of h, 5 and e) such that the 
number N(Ps,T) of eigenvalues of Ps in T satisfies 

|iv(Ps - r) -W^ vol(!, " (r),l - c ^ (m) 

with probability 

, C 

This is a restatement of Theorem 110.11 In Theorem I1U.3I we give a similar 
statement about the simultaneous Weyl asymptotics for all Ts in a family of sets 
that satisfy the assumptions of the above theorem uniformly. The lower bound on 
the probability becomes slightly worse but is still very close to 1 for suitable values 
ofe. 

The condition fll.lOp says that dT does not intersect the set of critical values of p 
and this is clearly not a serious restriction when T is contained in the interior of E. 
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However, we also would like to study the eigenvalue distribution near the boundary 
of S, and we then need a weaker assumption. 

Let r (s Q be open with C°° boundary. For z in a neighborhood of dT and 
< s, t <C 1, we put 

V z {t) = Vol {p E R 2 ™; \p{p) - z\ 2 < t}. (1.12) 

Our weak assumption, replacing fll.lOp is 

3k e]0, 1], such that V z {t) = 0(t K ), uniformly for z G neigh (<9r), < t < 1. 

(1.13) 

Here we have written in an informal way "neigh (dT)" for some neighbourhood of 
dT. Notice that ffTTOl) implies (TL~T5jl with k = 1. 

Generically, we will have {p, {p, p}} 7^ when p(p) 6 <9E and if we assume that 

at every point of p _1 (0), we have {p, p] 7^ or {p, {p,p}} 7^ 0, (1-14) 

then as shown in Example 112. 11 we have fl!.13[) with 5q = 3/4. (When p is analytic, 
we believe that (11.131) will always hold with some k > but have not consulted with 
experts in analytic geometry.) Under this more general assumption, we have 

Theorem 1.2 Assume $1.13\) and let 5 satisfy 

< 5 < h 3n+1/2 . 

Let N(P + SQu, T) be the number of eigenvalues of P + SQ^ in T. Then for every 
fixed K > and for < r <C 1 : 

\N(P + 5Q^ Y) - — ±— [ [ dxd£\ < (1-15) 
(27r/i) n JJp-i(r) 

(- + C K (r K + \n(-) [[ dxd£)), < r « 1, 



wzt/i probability 

> 1 _ 5 e -S(*r/0-» (1.16) 
r 

provided that 

/i K lni<e<l, (1.17) 


or equivalently, 

e'cf^ <5, C > 1, e < 1, 
implying that e > C/i K In r ; /or some C > 0. 
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This is a restatement of Theorem 1 1 2 . 3 1 and as explained after that theorem, when 
k > 1/2, the integral in the right hand side of (11.151) is 0(r 2K ~ l ) and it follows that 
we have Weyl asymptotics with probability close to 1, if we let r be a suitable power 
of e. To have the same conclusion when k < 1/2 we can assume that the integral is 
O(r ao ) for some ao > 0. 

Again we have a similar theorem about the simultaneous asymptotics for N(P + 
SQu, T) when T varies in a bounded family of domains satisfying all the assumptions 
uniformly. See Theorem II 2. 41 

The proofs follow the same general strategy as those of [7j with some essential 
differences: 

We do not use any non- vanishing assumption on the Poisson bracket - {p, p} . Instead 
we work systematically with the operators P*P and PP* and their eigenfunctions 
in order to set up a Grushin-problem. 

As in [7] we reduce ourselves to the study of a random holomorphic function, but in 
the present work this function appears as the determinant of the full operator (es- 
sentially) and we need to make some estimates for determinants of random matrices, 
and especially to prove that such a determinant is not too small with a probability 
close to 1. Those estimates were sufficiently elementary to be carried out by hand, 
but we think that future generalizations and improvements will require a careful 
study of the existing results on random determinants and possibly the derivation of 
new results in that direction. See the book [5] of V.I. Girko. 

Acknowledgements. We are grateful to F. Klopp for helping us to find some ref- 
erences. The first author was supported by a postdoctorial fellowship from Ecole 
Polytechnique. She also thanks Y. Colin de Verdiere for an interesting discussion 
around random functions. The second author is grateful to the Japan Society for 
the promotion of Sciences and to the Dept of Mathematics of Tokyo Unversity for 
offering excellent working conditions during the month of July, 2005. He also thanks 
E. Servat for a very interesting remark. We also thank the referee for many detailed 
remarks that have helped to improve the presentation. 

2 Determinants and Grushin problems 

Here we mainly follow [12] and give a more explicit formulation of one of the results 
there. Let Hi, TC2, Gi, G2 be complex Hilbert spaces and let Aj^ '■ Ti-k — ► Gj be 
bounded operators depending in a C 1 fashion on the real parameter t e]a, b[. We 
also assume that A^ k are of trace class and continuously dependent of t in the space 
of such operators. Here "over-dot" means derivative with respect to t. 
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Proposition 2.1 (flUtf) Assume in addition that A = (Aj^) : TCi x 7i 2 — > £1 x £/ 2 

/jas a bounded inverse B = {Bj^), and that A 2<2 and B\ x are invertible. (The 
invertibility of one of A 22 , -81,1 implies that of the other.) Then 

tr AB = tri 2 ,2^ 2 ,2 -trS^Sj,!. (2.1) 

Proof: We expand 

-^i- fc = ~~ ^2 ^j^A^B^, 

U fj, 

that are of the trace class too. In particular, 

— -Bi.i = ■Ei.i-Ai.i-Bi,! + 5i, 1^4^2-82,1 + -81,2^2,1-81,1 + -81,2^2,2-82,1- 

Rewrite the right hand side of (12. ip : 

tr i 2 ,2^2,2 ~ tr B ±}B 1A 
= tr A 2 ^ 2 A 2 l + tr A 1A B hl + tr A h2 B 2jl + tr B^B li2 A 2 ^B 1A + tr B^B h2 A 2 ^ 2 B 2A 
= tr A 2>2 A^ + tr A 1A B hl + tr A h2 B 2A + tr A 2) i-B li2 + tr B^B^ 2 A %2 B %1 
= tr A %2 {A~l + B 2:1 B^B li2 ) + tr A 1A B 1A + tr A 1>2 B 2A + tr A 2il B lj2 
= tr AB. 

Here we used the cyclicity of the trace and for the last equality the fact that 

-82,2 = ^2,2 + B 2t iB^lB li2 . (2.2) 

To check (12.21) we proceed by equivalences: 

(|2.2p -<=^ A 2 ^l = B 22 — B 2 ^B^lB 1<2 

-<==>- 1 = A 2j 2-B2, 2 — A 22 B 2 iB 1 1-81,2 
-<=>• 1 = (1 — ^2,1-81,2) + ^2,1-81,1-8^ Si, 2 

Here and in the following, we often denote the identity operator by 1 when the 
meaning is clear from the context. □ 

Consider the case TC\ = Q\ = TC, TC 2 = Q 2 = C N , 



A = V 



P i?„ 
R, 



S 



Assume also that P, V are invertible. (In the proposition we can permute the indices 
1 and 2 and think of P as A 2y 2-) We look for 

V = ( ~ ) : H x C N - H x C N , 

\R+ R+-J 

that is also invertible, i.e. we should be able to solve uniquely 

(2.3) 

i n + u -t- -ft + _'u_ = v + . 

Let 



Rewrite the first equation in (I2.3P as Pu = v — The general solution to that 

equation is 

u = E(v — RJu-) + E + v + , 
where v + , U- should be determined so that 

= E_(v - iLu-) + E_ + v + . (2.4) 

The second equation in ( 12.31) becomes 

v+ = R+E(v - RJa-) + R+E + v + + R + _u_. (2.5) 

Hence we get the following system that is equivalent to (12. 3ft : 



1 


Pu + 


- = v, 






R + u + R+- 


U- = 


v+. 


( E 


E + \ 


( p 


R\ 


U- 




\R + 


i 



E^ + v + - E-RJu- = -E-v, 

R + E + v + + (R + ^ - R + ER_)u^ = v + - R + Ev, 



[2.6) 



so (12.31) is well-posed iff 

TP TP TD \ 

: C 2N -> C 2N (2.7) 



E_ + ~ E -R- 

R+E+ Rj — R+ER- 



is invertible. 

Choose R + = tR + , R- = sR-, R + - = rid C N with s,t,r e C, and use that 
R+E + = 1, E_R_ = 1, R + E = 0, to see that the matrix (12.71) is equal to 



E_ + -s 
t r 



[2.8) 
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This matrix is invertible precisely when (s, t, r) belongs to the set 

{(s, t, 0); st ± 0} U {(s, t, r); r ^ 0, g <r(£_+)}. (2.9) 

Since P is invertible, we know that ^ <t(P |_). We can therefore find a C 1 -curve 

[0,1] 3 r h-> (s(T),t(r),r(T)) G the set (1231) . 

with 

(s(0),t(0),r(0)) = (1,1,0), (s(l),t(l),r(l)) = (0,0,1). 
This means that we have a C 1 deformation 

P M = {t { r)R + Vr)l') : « XC "^ XC " 
of bijective operators with 

V(0) = V, V(l) = J 

Applying (12. ip with the indices "1" and "2" permuted gives 

trVV- 1 = tr PP- 1 - tx EZ\{t)E- + {t) = -tr EZ\(t)E-+(t), 

where now "over-dot" means derivative with respect to r. If we integrate from r = 
to r = 1, we get with a suitable choice of branches for In: 

IndetP-lndetP = lndet £_ + (0). 

For this relation to make sense we also assume that 

P — 1 is of trace class. (2-10) 

Then for the original operator and its inverse we have 

In det P = lndet V + lndet P_+, (2.11) 

or equivalently, 

detP = detPdetP_+. (2.12) 
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3 General frame- work and reduction to trace class 
operators 

Let m > 1 be an order function on R 2n in the sense that there exist constants 
Co > 1, N > 0, such that 

m(p) < C (p - p) No m(p), Vp,p G R 2n , 

where we write (p) = a/1 + |p| 2 . We consider a symbol 

P(p; h) ~ + M/>) + ... in S(R 2 ", m), 

where 

S(R 2n ,m) = {a G C°°(R 2n ); |d^a(x,OI < C a m(x,£), V(x,£) G R 2n , a G N 2n }. 
Put 

£ = p(R 2n ), Soo = {z G C; 3pj G R 2n , j = 1, 2, 3, p^ -> oo, p(pj) ^ z, j ^ oo}. 

Assume 3z G C \ S, C > 0, such that 

- Zq\ > m{p)/C , WpeR 2n . (3.1) 

Then as pointed out in [7], for every z G C \ E, there exists C > such that 

b(p) - *| > m(p)/C, Vp G R 2 ", (3.2) 

and for every z G C \ Eoo, there exists C > such that 

\p{p) -z\> m(p)/C, Vp G R 2n with |p| > C. (3.3) 

Let fl d= C\S 0O be open simply connected containing at least one point zq G C\S. 

Lemma 3.1 For every compact set K G Q, there exists a smooth map k : Q\{zq} — > 
Q\{zo} with k(z) = z for all z in a neighborhood ofdQ, such that K,(Y,r\Q)r\K = 0. 

Proof: Q is diffeomorphic to the open unit disc D(0,1) in such a way that Zq 
corresponds to 0. Now consider k : D(0, 1) \ {0} — > D(0, 1) \ {0} defined by k(z) = 
f(\z\)z/\z\, where / is a smooth function on ]0, 1] with 1 — e < f(r) < 1, such that 
f(r) = r for 1 — r < e/2. Choosing e > small enough and conjugating with the 
diffeomorphism above we get the desired map k. □ 
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Let n <e f2 be open. Take k as in the lemma with K containing the closure of 
Q. Extend k to be the identity in C \ Q and put p = nop. Then p(p) — z is elliptic 
in the sense of (13. 2p . uniformly for z G f2 and 

p-peC™(R 2n ). (3.4) 

Put 

P = p + hpi + h 2 p 2 + ... . 

Now pass to operators and denote by the same letters symbols and their h-Weyl 
quantizations. We shall consider P as a closed operator: L 2 — > L 2 with domain 
H(m) := (P — Zq)~ 1 L 2 (see [7]). From the discussion above, we get 

• For every compact set K C C \ S, we have cr(P) fl K = 0, when /i > is small 
enough. 

• cr(P) C\Q is discrete when h > is small enough. 

• cr(P) fl Q = when /i > is small enough. 

In view of the last property and (13. 4ft . we also have (for h > small enough), 

Proposition 3.2 For z EVt, we have that 

P{z) := (P - z)-\P -z) = l + K(z), 

where K(z) is a trace class operator. Moreover, 

z E a(P) & E cr(P(z)). 

Notice that K(z) = (P — z)~ 1 {P — P) is the quantization of a symbol belonging 
to the intersection of S(m) for all order functions m. 

4 Some functional calculus 

Let P = 1 + K, K E Op h (S(m)), where m E C°°(R 2n ; ]0, oof) is an integrable 
order function. We also assume that K = k + hki + ... in S(m) on the symbol 
level. We shall review some functional calculus for Q = P*P and more generally 
for a selfadjoint operator Q > with Q ~ q + hqi + ... on the symbol level, with 
Q- 1 G S(m), q>0. 

Let -0 £ Co°(It)- For /i < « <C 1 we shall study the properties of ij)(ar l Q). 
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To this end we shall consider a~ l Q as a symbol with /i/aasa new semiclassical 
parameter, after a suitable dilation in phase space. More precisely, we make the 
change of variables 

x = a^x, D x = a~?Dx 

and write 

-Q(x, hD x ; h) = -Q(a*(x, -D*)\ h), (4.1) 
a a a 

with symbol a~ 1 Q(a 1 ^ 2 (x, £); h) for the /i/a-quantization. The lower order terms 
are O{h/o) uniformly with all their derivatives, so we shall just look at the leading 
symbol 

a 

where we dropped the tildes on the new variables. The (new) associated order 
function will be 

m(x,Q :=l+ q ^ M) >l. (4.3) 

a 

We have 

Vm = (V g )(a*(s,fl) < c qHaHx,Z)) < Cm(x,OK 
V 2 m = 0(1), 

so by Taylor's formula, 

m(p) = m(/i) + 0(l)m((i)*\p - fi\ + 0(l)\p - p\ 2 , 
and since m{p) > 1: 

m(p)<C(p-/i) 2 m(/i). (4.4) 

Hence m is an order function, uniformly with respect to a. 
Similarly, we have the improved symbol estimates 

V^ = 0(l)m*, (4.5) 

v2 ^!p) = (46) 

a 

V fc ^M = fc > 2 . (4.7) 
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In particular, we have the standard symbol estimates 

v fe^P) = o(l)m(p). (4.8) 
a 

It is therefore clear that we can apply the functional calculus in the version of [9J 
(see also [2]), to see that if ip £ Cq°(R), and if we interpret Q/a as the right hand 
side of ( 147TT) . then 

ii>(a- 1 Q)=Oj>H s (f), (4.9) 

where 

oo , 

/ = E(-r/^,f), in^m" 1 ), (4.10) 

V a 

with/ (^,0 = V'(a- 1 g(« 1/2 (J,0)), 

= E a^^^a)^ ^"" 1 ^" 172 ^^))), aj> e S(l). (4.11) 



Proposition 4.1 Let m = m a (x,£) be an order function, uniformly with respect to 
a, such that m(x,£) = 1 for a~ 1 q(a 1 ^ 2 {x,^)) < supsupp-?/; + l/C , for some C > 
that is independent of a. Then holds in S(fh), for h and h/a sufficiently 

small. 

Proof: Write q a = g(a 1 / 2 (x, £))/ a > Q a = a~ 1 Q(a 1 / 2 (5;, -D^); h) and drop the tildes. 
Let q a G S{m) be such that supsuppV' + l/(5C) < q a , and be equal to q a when q a > 
supsupp?/> + 2/(5(7). Let Xa G 5(1) be equal to 1 when g Q < supsupp-^ + 3/(5(7) 
and equal to when g a > supsupp ip + 4/(5(7). We use the same symbols q a , q a , 
Xa to denote the h/a quantizations. 

Let ip be an almost holomorphic extension of ip and recall the Cauchy-Green- 
Riemann- Stokes formula, in the operator sense ([9], [2J, [3]): 

^(<? Q ) = ^ J(z-q a y ll ^L(dz) } 

where L(dz) denotes Lebesgue-measure. For z in a neighborhood of supp ?/>, we write 
(z-q a y l = (z-g a ) _1 ox a + (2;-g a ) _1 o(l-x Q )-(2;-g a ) _1 (g Q -g a )(2;-g a ) _1 (l-Xa). 
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Then, since the middle term is holomorphic near the support of ip, 

1 f d%f> 
4>(q a ) = 4>(q a ) °Xa~ - (z- qa)' 1 ^ - q a )(z - g a ) _1 (l - Xa)^=L{dz). 

Here the symbol of ip(q a )ox a has the asymptotic expansion ( 14.10(1 in S(fh), thanks to 
the extra factor \ a and with the same terms given by (14. lip . For z G neigh (supp ip), 
(z—q a )~ l G Op (1/m) depends holomorphically on z and thanks to the factor q a —q a , 
whose support on the symbol level is separated from that of 1 — Xa by some fixed 
positive distance, we know that (q a — q a )(z — g Q ,)" 1 (l — Xa) G (h/a) N Op (S(fh)) for 
any N > and any fh as in the proposition. Combining this with the estimates for 
the symbol of (z — qa)" 1 from the Beals lemma as in [9] (see also [2], Proposition 
8.6), we get the proposition. □ 

We next apply the functional result to the study of certain determinants. Let 
X G Q° ([0, +oo[; [0, +oo[), x(0) > and extend x to C£°(R; C) in such a way that 
x(t) > near and t + x(t) ^ 0, Vt G R. We want to study lndet(Q + ax(a _1 Q)), 
when /i C a < 1. Let us first recall from [TU] that if Q = Op h (q) with q G S(l), 
q > 0, q — 1 G S(fh), where fh is an integrable order function, then 



lndetQ = — —^{JJ \nq(x,t)dxd!; + 0(h)). (4.12) 

In fact, let Q t — (1 — t)l + tQ, so that Qo = 1, Qi = Q. Then by standard elliptic 
calculus, with q t = (1 — t) + tq, we have 

-^lndetQt = txQ^—Qt 
dt dt 



(2nh) n 



[JJVj t QtdxdC + 0(h)) 



1 ^ I I \nq t (x,£)dxd£ + 0(h)), 



(2nh) ny dt 

and integrating from t = to t = 1, we get (14.121) . 

For a = ai > fixed with «i <C 1, this applies to Q + aix(oii 1 Q) and we get 

lndet(Q + a lX (a^Q)) = J^j^i j j ln (<? + a lX (a^q))dxd^ + 0(h)). (4.13) 
We have for t > and E > 0: 
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with 

1 ' E + x(E) 
Now for h -C a < £ < ai, we get from Proposition 14. II by dilatation: 

-|mdet(Q + £x(r 1 Q)) =trrV(* _1 Q) (4-14) 

J! (i + dist ((x, g; supp x(^P-)T N dxdtj , 

where < x G C^°(R) is equal to one on some interval containing [0, sup supp if}] 
and the last term is coming from the "remainder" in the asymptotic development 
P~ID|) . 

We are interested in the integral of this quantity from £ — a to £ — ot\. Let us 
first treat the leading term 

t „ is \ i>( ^Mi )dSdi = i // ^(^)drfe 



v 2vr/i y JJ t t ' {2nh) n J J £ v £ 

and integrating this from £ = a to £ = cti, we get 

1 1 ' \n{q + t X {-))dxdt] ai . (4.15) 



.{2nh) n JJ y H 
The second term in (14. 14ft is 



t=a 



r 2 

Integrating this from £ = a to £ = a±, we get 



< 0(1)^ // r/,7/4. 



g(z,£)<C 2 t 



^(l)/^™ /TV \dtdxdi (4.16) 

0(1)^-" / / [ ? \ — -r - — I 



9(*.0<Cai \max (a,q(x,£)/C) Qfi 
7 \" ( 



< 0(l)/i~ n / / h - dxd£ 
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When estimating the third term in (14.141) we consider separately the regions 
1(^)01 — C an d 1(^)01 > C f° r some large C > 1. Consider first the region 
< C. Put 

dt(5,f) = dist ((£,£), {(y,7/); -q{t^(y,r))) < 6}), C = supsuppx- 

For (y, rj) with g(£2 (y, rf))/t< C, we have V(g(^ (y, 77))/^ = 0(1) and V 2 (g(t^ (y, 77))/^ 
0(1), so by Taylor expanding at (y,rj) we get 

\q{£ (x, 0) < o(i)(i + d t (x, g + o 2 ) < o(i)(i + 0) 2 . 

The contribution to the third term in (I4.14p from t 1 ^ 2 \ (x, £) | < C is therefore 

N {l){^r\l! {l + \q{£{x^))Y N dxdl (4.17) 



= ^(1)^)^ // ( 1 + jq(x,t))- N dxd£. 
We integrate this from t — a to a±, so we want to estimate 

/•ai -I 
,M ' 1 



h M / : dt 

If q < a, we get 

If a < q < oti, we get 



t M+iy t > t M+iy t > j a t M+i qN y q > > 

with the symbol ^ indicating same order of magnitude. Choose N = M + 1, to get 

-(-) . 

For «x < g, we get 
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with the same choice of N. Thus the expression (14.171) is 

c»(i) rr / h ^ M 



dxd£, VM > 0. 



We next look at the contribution to the last term in (14.141) from the region 
tz\(x, £)| > C, which is 

°^ilTl I A (l + \(x,0\)~ N dxdi= 0(l)(Vt* VM,iV, 

= 0(l)/i M , VM. 

Summing up, we have proved 

Proposition 4.2 Let x e C^(R) uritfi x(0) > and x(t) > for t > 0. T/jen /or 
/i < a < 1: 

In det(Q + ax{a' l Q)) = (4.18) 

[ [ ln(q + ax(^))dxd^ + 0(1) [[ ) dxdi)+0{h™). 

{2Txh) n \JJ a JJ\{x,o\<ca + q{x^) ) 

Most of the proof was based on (14.141) , of which the second part is valid for any 
if) G (7q°(R). The estimates leading to the preceding proposition, also give 

Proposition 4.3 Let x £ Cq?(R), and choose x G C^(R; [0, 1]) equal to 1 on an 
interval containing and supp x- Then for < /i < a < 1, we have 



tr X (a- l Q) = J^ijj x( q -^)dxd^ + OQ Jj x( q ~^)dxd^ (4.19) 

+ ^,A/(1)(-) M // (l + ^)- Ar rfx^ + 0(^° O )). 



Put 



-) M (f (l + ^- N - 

V(t) = [[ dxdti, < t < \, (4.20) 

so that < V(t) is an increasing function. By introducing an assumption on V(t), 
we shall make (I4.19P more explicit and replace (I4.18P by a more explicit estimate. 
We assume that there exists k e]0, 1] such that 

V(t) = 0(l)t K , < t < -. (4.21) 
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The first integral in (I4.18P can be written 

J I ln(g + ax(~))dxd£ = J J \n(q)dxd^ + 0(1) J J In ~dxd£ 



q<Cc 
Ca 



j j \n(q)dxd£ + 0(1) J In ^dV(q) 



so (14.181) gives 



\ndet(Q + a X (-Q)) (4.22) 
a 

i , rr _ r Ca . i „, , h 



} n ([[Hq)dxdZ + 0(l) [ a \n-dV(q) + 0(l) f* -^-dV(q) + 0(h°°)). 
[2'Kh) n JJ Jo q io a + 9 

Similarly (14.191) can be written 

trxia-'Q) = j^j-t [ X (-)dV(q) + 0(-) [ X (-)dV(q) (4.23) 
{27ih) n V J a a J a 

+O n ,m{1){-) M Hi! + -r N dV(q) + 0(h)). 

a Jo a J 

In particular, the number N(a) of eigenvalues of Q in [0, a] satisfies 

N(a) < 0(l)(h- n f \l + ^-)- N dV(q) + h l - n ). (4.24) 
Jo a 



Proposition 4.4 Under the assumption fiJJHty , we have 



pa. 

/ hx{q)dV{q) = 0(a K \na), (4.25) 
Jo 

ai , (0(a"±), for K <l, 

—dV(q) = (4.26) 
o a + Q 

yO(hlia-), when k = 1, 



N(a) = 0(a K h~ n + h 1 - 71 ). (4.27) 



19 



Proof: This follows by straight forward calculations, starting with an integration 
by parts. 

f a in(q)dV(q) = [Hq)V(q)]% - f -V(q)dq = 0(a K \n(a)), 
Jo Jo 1 

/ -^-dV(q) = l—-V(q)}^ + / — - 2 V( q )dq 
Jo a + q a + q J [a + q) 2 

f ai ha K 
= 0(h) + 0(1) / j^-^dq 
Jo (« + Q) 

pai/a 

= 0(h) + OiTjha"- 1 



Q 



(1 + g) 2 * 

(D^ha*- 1 ), 0<k<1, 

0(Mn±), k = 1. 
To get P?| , we use ffl~24l) and the estimate 

r\i + ^r N dv(q) = i(i +^r N v(q)]?+N / Q1 (i + -)- (iV+1 ^(?)- 

Jo a a Jo a a 

POO 

= 0(1)0* + 0(1) / (I +q)-( N+l) q K dqa K 
Jo 

= 0(l)a K . 



□ 



Since we will always assume that /i < a < 1, and since k < 1, we can simplify 
g]23) to 

N{pt) = 0(l)a K h~ n . (4.28) 

Proposition 4.5 Assume \J^.21^ . Under the assumptions of Proposition we 
have 

lndet(Q + ax(a _1 Q)) = , \ , ( 1 1 Mq)dxd£ + 0(l)a K In a). (4.29) 

(2-Kh) n J J 

Under the assumptions of Proposition we have 

txx^Q) = -\-( [ [ X (^ L )dxdi + G(l)a K -). (4.30) 
For < ft < a < 1, the number N(a) of eigenvalues of Q in [0, a] satisfies 
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Notice that when Q > 0: 

IndetQ < lndet(Q + 

so (I4.29P with a = Ch, C 3> 1, gives an upper bound which is more precise than 
the one in [TO] . 

5 Grushin problem for the unperturbed operator 

In Section [3] we introduced the operator 

P(z) = l+K{z), K(z)eO Vh (S{m)) 

where m is an integrable order function, so that K(z) is a trace class operator. 
Here P(z) depends holomorphically on 2 6 O d C, where Q is open. Also recall 
that P(z) = (P — z)~ l (P — z). We are interested in the spectrum of small random 
perturbations of P; P$ = P + SQ. Correspondingly, we get Pg(z) = (P — z)~ 1 (Ps — 
z) = 1+Ks(z), and the main work in later sections will be to study | det(l+ Ks(z))\. 
The upper bounds will be fairly simple to get, and the delicate point will be to get 
lower bounds. As a preparation for this more delicate step, we here study a Grushin 
problem for the unperturbed operator P(z). In this section z G Q will be fixed and 
we simply write P instead P(z). 

Let ei,e2, ... be an orthonormal (ON) basis of eigenvectors of P*P and let < 
Ai < A2 < ... be the corresponding eigenvalues. (Strictly speaking, if we want the 
eigenvalues to form an increasing sequence, the set of indices j should be of the form 
J = Jo U Ji U J2, with 

• Jo = N or a finite set, < Xj < 1, for j G Jq, 

• J\ = N or a finite set, Xj — 1, j G J\ 

• J 2 = — N or a finite set, Xj > 1, j G Ji- 

We will only be concerned with finitely many indices from Jq.) 

Since P = P(z) is a Fredholm operator of index zero by Proposition I3.2[ we 
know that PP* and P*P have the same number iVo of eigenvalues equal to 0. Let 
fx, /tv be an ON basis of ker (PP*). For j > Nq, we have Xj > and Pej is an 
eigenvector of PP* with eigenvalue A^: 

PP*P ej = XjPtj. 
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Using standard notations for norms and scalar products, we put fj = \\Pej\\ 1 Pej. 
Then {fj}j^j is an ON system of eigenvectors of PP*, with PP*fj = Xjfj. Let 
/ G L 2 (R n ) with (f\fj) = for all j G J. Then (P*f\ej) = U\ Pe j)- If 3 < N o, we 
get (P*f\e 3 ) = (f\0) = 0, and if j > N + 1, we get (P*/|e,) = \\PejUf\fj) = 0. 
Hence P*f = 0, so / G ker(PP*) and hence / = is zero since ker(PP*) is the span 
of /i, /jv - We conclude that is an ON fraszs of eigenvectors of PP*. By 

construction, we have Pej = Wjfj with < Wj = ||Pej||. Then 

w] = (P*Pe j \e j ) = \ j , 

so Wj = y/Xj and it follows that, 

Pe, = V~,fr (5-1) 

P*fj = v'V, (5-2) 

for all j £ J. 

Let < a < 1 and let iV = N(a) be given by 

\ j < a <=> j < N(a). (5.3) 



Define 
by 

and put 



R + : L 2 (R n ) -> C N , P_ : C w -> L 2 (R n ) 

AT 

P+w = y/a((u\ej))* =1 , R-U- = Va ^ u -0)fji 



P P_ 
R, 



: L 2 (R n ) xC^ L 2 (R n ) x C*. (5.4) 



N 



v ( u \ = (Y.j V\jUjfj + Ei Vau-U)ff\ 

W V (^)f=i J ' 

and we conclude that 

f - / \ . . F7\l\ -m w , 

\i) (5.5) 

N<jeJ 

N 



detPi = <ni**(^ n 

i v v 7 N<jeJ 

n 

N<jeJ 

= or* Y\ max(\/a, y/Xj). 
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Notice that 

| det P| = Yl (5-6) 
j 

Let Sj(k) = 5j t k, 1 _ j,k < N. Then V maps Cej x C<5j to Cfj x C£j and has 
the corresponding matrix 



o 4= 



A, ) 



The inverse is given by 



i 

so if v — ^2jVjfj, v + = J2 1 v + (j)8j, then (writing as before £ for P -1 ) 
c ( v \ 75^ + * E " V+U)ei \ = (E E + \( V] 

' ^ ~~ { Ef frA - Ef ' ~ ^" £ ^ ^ ' * 



where 



1 * 

P+^+ = ~f= ^u+0*) e i» ( 5 - 8 ) 



1 ? 

P-+ = -^diag(A/A~), 

||p||jp + ||,||i5_||, IIJ5_j.II < 



1 



a 



From fl5.5l) - fl5.8p . we see that 

| det P | = |detP||detP_+|, (5.9) 

as we already know from (I2.12p . 

We next study | detP|, when ZiCaCl. The formula (15.51) can be written 

|detP| 2 = a 7V detl Q (P*P), (5.10) 

where l a (t) = max(a,t). Let x £ Co°([0, 2[; [0, 1]) be equal to 1 on [0, 1]. Then for 
t > 0, 

t + jx(-)<Ut)<t + a X (-). (5.11) 
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In the following, we assume that Q = P*P satisfies the assumptions of Section 
HI including (14.211) . and choose h <C a <C 1. Then we know that 

N{a) = 0{a K h~ n ), (5.12) 

and Proposition 14.51 in combination with f l5.10l) -( 15.12l) show that 

ln|detP| 2 = 1 ( JJ\n(q)dxd£ + 0(l)a K \na). (5.13) 

As noticed after Proposition 14. 5[ we also have the upper bound 

lndetP*P< . \ . ( I [ ln(q)dxd£+ C(l)a K ln-). (5.14) 
{2nh) n J J a 

6 The Hilbert-Schmidt norm of a Gaussian ran- 
dom matrix. 

Let a(u) be a complex Gaussian random variable with density 



1 



7T<7 2 



L(da), L(da) = dRea dlma, (6.1) 



that is a A/"(0, <7 2 )-law with a 2 denoting the variance. The distribution of |a(u;)| 2 is 

^da = -e- r/s H(r)dr, (6.2) 
s 

where s = a 2 and H{r) denotes the standard Heaviside function. Notice that 

\\\a\ 2 \\ L , = (\a\ 2 ) = a 2 . 

Let <x,(c<j), j = 1,2, ... be independent random variables as above with variance 
crj and assume for simplicity that o\ > <jj for all j. We also assume that 

oo 

$>!<oo, (6.3) 
i 

implying the a.s. convergence of l^'^)! 2 - 
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We want to estimate the probability that |a!j(u;)| 2 > a. The probability 
distribution of l^jO^)! 2 i s ec L ua l to (/ii * /i 2 * ....)dx, where fij is given in (16. 2|) 



with s = Sj = Uj, so that 



i 



The Fourier transform of fij is given by 

Mp) = irk? (6 - 5) 

which has a simple pole at p — i/sj. The probability that we are after, is 
f°° 1 f 00 — 

/ (/ii * /i 2 * ...)dr = — / YlCflj(p))l\^[(p)dp, (6.6) 

•/a 1 

by Parseval's identity. Here 

so the probability (16.51) becomes 



2 



The assumption (16.41) implies that the infinite product converges away from the 
poles i/sj. For p in a half plane Im p < 6 < we have 



1 

< 



1 + ((l-b Sl ) 2 + sl(Re P y 



2 ' Jr 2 J 2 

where Co is a universal constant appearing in the estimate, 

< e Co \ 0<t< 



l-t~ ~ ~ 2 
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Shifting the contour in (16. 8p from R to R + ib and choosing b = l/(2si), we can 
estimate the probability (16. 6p from above by 

C °° 1 
C{s x )exp[^Y^s j - — a], (6.9) 

where Co > is the universal constant introduced above and C(s\) can be chosen 
uniformly bounded on any compact subset of ]0, +00]. 

Remark 6.1 W. Bordeaux Montrieux has used a more elementary argument in 
the case of real matrices, by means of the Markov-Chebyschev inequality. With P 
denoting the probability and () the expectation values, it gives for every a > 0: 



OO 



a p(j2 K-nr >*)<(£ i^-hi 2 ) = E<i a »i 2 > = E4 ( 6 - 10 ) 
1 111 

We will prefer (16. 9p however, since it gives an exponential decay with respect to a. 

Remark 6.2 If Q = (aj t k(^))j,k£'N is a random matrix where aij^uj) are indepen- 
dent jV(0, a^ k ) laws, and 

J>f, fe <oo, (6.11) 

then (16. 9p gives an estimate on the probability that the Hilbert-Schmidt norm is 

>a 1 / 2 : 

P(||K,(o;))||^>a)<C( Sl )ex P [^ £ aj fc - i-a] (6.12) 
where C , C(si) is the same constants as in (I6.9P and s\ = maxer? fc . 

7 Estimates on determinants of Gaussian random 
matrices 

Consider first a random vector 

t u(co) = (a 1 (u),...,a N (u))EC N , (7.1) 

where ati, are independent complex Gaussian random variables with a A/"(0, 1) 
law and to is the random parameter living in a probability space with probability P. 
The law of a>j, i.e. the direct image of P under otj, is given by 

K-)*(P) = -e"l z l 2 L(^) =: f(z)L(dz) (7.2) 

7T 
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and L(dz) = L c [dz) is the Lebesgue measure on C. 
The distribution of u is 

u,(P) = \e-\ u \ 2 L cN (du). (7.3) 

If U : — > is unitary, then C/w has the same distribution as u. 

We next compute the distribution of |it(o;)| 2 . The distribution of |<x,'(c<j)| 2 is 
fi(r)dr, where 

fjtM = -H(r)^-e- r = e- r H(r), 
dr 

where H(r) = l [0:Oo[ (r). We have p(p) = 

We have |w(u;)| 2 = \aj(u)\ 2 and since la^u;)! 2 are independent and identi- 
cally distributed, the distribution of |w(u;)| 2 is fi*...*fidr = fi* N dr, where * indicates 
convolution. For r > 0, we get by straight forward calculation the x%n distribution 
(for the variable 2r) 

N— 1 — r 

^ N dr = J^yH(r)dr. (7.4) 

Recall here that 

oo 

JV-1 -r 



r N - x e~ r dr = T(N) = (N-l)\, 

so [i* N is indeed normalized. 

The expectation value of each |aj(to>)| 2 is 1 so: 

(\u(u)\ 2 ) = N. (7.5) 

We next estimate the probability that |u(u;)| 2 is very large in a fashion that is 
slightly different from that of Section [61 It will be convenient to pass to the variable 
ln(|tt(u;)| 2 ), which has the distribution obtained from (17. 4p by replacing r by t = lnr, 
so that r = e*, dr/r = dt. Thus ln(|-u(u;)| 2 ) has the distribution 

:H{r)— = — — —rdt =: u N (t)dt. (7.6) 



(iV-1)! v ' r (N — 1)! 
Now consider a random matrix 

(ui...u N ) (7.7) 
where Uk(w) are random vectors in (here viewed as column vectors) of the form 
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and all the aj^ are independent with the same law (17.21) . 
Then 

det(iti u 2 ...un) = det(«i u 2 ...un), (7.8) 

where Uj are obtained in the following way (assuming the Uj to be linearly indepen- 
dent, as they are almost surely): u 2 is the orthogonal projection of u 2 in the orthogo- 
nal complement (wi) -1 , M3 is the orthogonal projection of u% in (ui,^)" 1 = (ux,u 2 ) 1 ', 
etc. 

If Ui is fixed, then u 2 can be viewed as a random vector in C^" 1 of the type 
(17.11) . (17. 2p . and with Ui,u 2 fixed, we can view u% as a random vector of the same 
type in C^" 2 etc. On the other hand 

I det(u! u 2 ...un)\ 2 = \ui\ 2 \u 2 \ 2 ■ .. ■ \u N \ 2 . (7.9) 

The squared lengths |wi| 2 , \u 2 \ 2 , |mat| 2 are independent random variables with dis- 
tributions fx* N dr, ix*^ N ~^dr, fidr. This reduction plays an important role in [5|. 
The following lemma will not be used directly. 

Lemma 7.1 Let a, (3 > be independent random variables with distributions // a (r)— ; 
Pp( r )^r- Then the product ct(3 has the distribution /i a| g— , with 

Ha/3 = V4pp ■= M~ l {{Mp, a ){M^)). (7.10) 

Here 

dr 



Mfx(r) = J r- iT fi(r) 

is the Mellin transform of \i. 



r 



Proof: Recall that the Mellin transform of \x{r) is the Fourier transform of /x(e*); 
r = e*, r~ x dr = dt. The distribution of In a is related to that of a by the same 
change of variables p a (r) — — > n a (e l )dt = u a {t)dt. Since multiplication on the 
Fourier transform side corresponds to convolution, (I7.10p is equivalent to the fact 
that the distribution of the sum of two independent random variables is equal to 
the convolution of the distributions of the two variables. □ 

The proof also shows that the multiplicative convolution in the lemma is given 

by 



Pa%Pp(r) 



[°° Ha(-)^(P)-. (7.11) 
Jo P P 



As already mentioned we shall not use the lemma directly but rather its proof by 
taking logarithms and use that the distribution of the random variable In | det(«i u 2 ...un 
is equal to 

[y x * v 2 * ... * VN)dt, (7-12) 
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with Vj defined in (17. 6p . 
We have 

u N (t) < u N (t) := - J— e ™ 
Choose x(N) G R such that 



u N {t)dt = I. (7.13) 



More explicitely, we have 



1^*^ = 1, x(iV) = lln(iV!) = llnr(iV+l). (7.14) 



Using Stirling's formula, 



(N-iy._r { N ) _ e _ NNN _ ki+ , 



v 7 ^ ' N 

we get 

<N) = I(Iln(27r)-(iV + l) + (iV + i)ln(iV + l) + 0(i)) (7.15) 
= ^{(N + l)\nN-N + C + O(±)) 

where C = (ln2?r)/2 > 0. 

With this choice of x(N), we put 

PN(t) = l]-oo,x(N)](t)u N (t) , 

so that px(t)dt is a probability measure "obtained from v^{t)dt, by transfering mass 
to the left" in the sense that 



fu N dt < J fp N dt, (7.16) 

whenever / is a bounded decreasing function. Equivalently, 

9 * < 9 * Pn, 
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whenever g is a bounded increasing function. Now, for such a g, both g * vn and 
g * Pn are bounded increasing functions, so by induction, we get 

g * v\ * ... * i>n < 9 * Pi * ••• * Pn- 
In particular, by taking g = H, we get 

(z/x * ... * v N )(s)ds < / (pi * ... * p N )(s)ds, t e R. (7-17) 



We have by f[T~i4j) 

p}v(r) = r^ ,^ 1 — e^-^rft = — 1 , e M»(*)-fcc(*)r ( 7,i 8) 



, x (JV-1)! ' (iV-l)l(iV-zr) 



1-^ 



This function has a pole at r = —iN. 
Similarly, 



W](r) = — ^ e " aT - ( 7 - 19 ) 



By Parseval's formula, we get 

/a -y /*oo 

pi * .. * p N dt = — / «F(p! * .. * p 7V )(r)J r l ] _ 00ia ](r)dt 

oo ^ 7I " ./— oo 



(7.20) 



1 /-+O0 • N -i 

- / e -^-y)-«)^!_TT_J_d T . (7.21) 
2W-00 r-i0Al(l-^) 1 J 



We deform the contour to Imr = —1/2 (half-way between R and the first pole in 
the lower half-plane). For j > 2, we use the estimate 

when Imr = —1/2. Hence, 

N , .V 



filial <«p(5E(i + ^))<o^- 



3 2 
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It follows that for a < ^ x x(j) : 

r a i N 

/ Pl * .. * PN dt < CNhxp - «))■ I 7 - 22 ) 

In view of (I7.17p . (I7.20p the right hand side is an upper bound for the probability 
that In | det(ui...U]y)\ 2 < a. 

From the formula ( 17. 151) . we get for some constants C\, C 2 G R: 



N 



> Ci + iN + ^lnN -2N + -(\nN) 2 + Co\nN > C 2 + (N + -)\nN -2N. 

1 

(7.23) 

Hence, for a < C 2 + (N + §) IniV - 2iV, 

P(ln I det(Mi...M7v)| 2 < a) (7.24) 
< CN*exp[~(C 2 + (N+-)kiN-2N-a)] 

= Cexp[-^(C 2 + (N-^)lnN-2N-a)}. 

We shall next extend our bounds on the probability for the determinant to be 
small, to determinants of the form 

det(L> + Q) 

where Q = {u\...un) is as before, and D = [dx...d^) is a fixed complex N x N 
matrix. As before, we can write 

I det((di + ui)...(dN + u N ))\ 2 = \di + u\\ 2 \d 2 + u 2 \ 2 ■ .. • \d N + u N \ 2 , 

where d 2 = d 2 (ui), u 2 = u 2 {ui,u 2 ) are the orthogonal projections of d 2 , u 2 on 
id\ + ui)- 1 , <i 3 = ds(ui,u 2 ), u 3 = u 3 (ui,u 2 ,U3) are the orthogonal projections of d 2 , 
u 2 on (di + Mi, d 2 + u^ 1 - and so on. 

Let i>^ N \t)dt be the probability distribution of In \d + u\ 2 , when d G is fixed 
and u G C N is random as in flUD, O- Notice that u^ N \t) = v {N \t) is the density 
we have alreay studied. 

Lemma 7.2 For every a G R, we have 

f vf\t)dt< f u^ N \t)dt. 

J— 00 J— 00 
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Proof: Equivalently, we have to show that P(|d+ u\ 2 < a) < P(M 2 < o) f° r every 
a > 0. For this, we may assume that d = (c, 0, 0), c > 0. We then only have to 
prove that 

P(|c + Reui| 2 < b 2 ) < P(|Re«i| 2 < b 2 ), b > 0, 
and here we may replace P by the corresponding probability density 

1 _ 2 

fi(t)dt = —^e~ l dt 
for Re/ii. Thus, we have to show that 

1 1 e -t 2 dt < -L f e~ t2 dt. (7.25) 



/7r J\c+t\<b V 71 " J\t\<b 

Fix b and rewrite the left hand side as 



The derivative satisfies (recall that c > 0) 



e~ t2 dt. 



r( c ) = -L( e -( b+c f - e^ c > 2 ) < 0. 



'IT 

hence c i— > 7(c) is decreasing and ( 17.25ft follows, since it is trivially fulfilled when 
c = 0. □ 

Now consider the probability that In | det(D + Q)\ 2 < a. If Xa(t) = H(a — t), 
this probability becomes 

P(dui)...P(duN) x 

Xa(hi \di + Mi| 2 + In |d 2 (wi) + M2(mi,M2)| 2 + ••• + In |djv(ui, Wjv-i) + %(mi, .., wat)| 2 ). 
Here we first carry out the integration with respect to u^, noticing that with the 
other Ui,..,«jv-i fixed, we may consider djy(ui, .., Mjv-i) as a fixed vector in C ~ 
(c?i + Mi, c?iv_i + mjv^!)- 1 and % as a random vector in C. Using also the lemma, 
we get 

P(ln|det( J D + Q)| 2 < a) 



v \ (tN)dtNP(duN-i)--P(du 1 ) x 



< 



Xa(ln |di + Mi| 2 + .. + In \d N -i(ui, ..,u N ^ 2 ) + «at-i(«i, ..,mat_i)| 2 + t N ) 
J .. J v {1 \t N )dt N P{du N ^i)..P{dui) x 

Xa(ln \di + wi| 2 + .. + In \d N -i(ui, ..,u N - 2 ) + u N -i(ui, ..,mat_i)| 2 + t N ). 
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We next estimate the u^-i- integral in the same way and so on. Eventually, we get 
Proposition 7.3 Under the assumptions above, 

P(ln|det( J D + Q)| 2 <a) < J .. J Xa ( tl + ... + t N )^\t N )v^ {t^)..^^) 

= P(ln|detQ| 2 < a). 

In particular the estimate \1-2J$ extends to random perturbations of constant ma- 
trices: 

P(ln | det(£> + Q)\ 2 < a) < Cexp [~(C 2 + - 1) InN - 2N - a)}, (7.26) 

_ _ 

when a<C 2 + (N + l)\nN- 2N. 

8 Grushin problem for the perturbed operator 

Let P be as in Section [3j Let < m, m < 1 be square integrable order functions on 
R 2n such that fh or m is integrable, and let S G S(m), S G S(fh) be elliptic symbols. 
We use the same symbols to denote the /i-Weyl quantizations. The operators S, S 
will be Hilbert-Schmidt with 

||£||hs, ||3||hs - h~%. 

Let ei, 62, and ei,e"2, ... be orthonormal bases for L 2 (R n ). Our random perturba- 
tion will be 

= S ° a oA u )^k ° (8.1) 

where at^ are independent complex 7V(0, 1) random variables. See the appendix, 
Section [13] for a general discussion. 
Consider the polar decompositions 

S = DU, S = UD, (8.2) 

where U, U are unitary pseudodifferential operators with symbol in S(l) and D, 
D are positive selfadjoint elliptic pseudodifferential operators with symbol in S(m) 
and S{m) respectively. After replacing ej by Uej and ej. by U*ek, we get with the 
new orthonormal bases that 

Qco = D aj^i^ejel o D, (8.3) 
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Now as we recall in the appendix (Section [15]) . we may replace the bases ej and ej by 
any new orthonormal bases we like, if we replace the aj tk (cu) by a new set of random 
variables (that we also denote by aj±) having identical properties. If we choose ej 
to be an orthonormal basis of eigenfunctions of D and similarly for e}, then we get 

Qu; = ^2s j s k a jjk (uj)e j e f k , (8.4) 

where "Sj > and Sj > are the eigenvalues of D and D respectively, i.e. the 
singular values of S and S. 

We are then precisely in the situation of Section noting that s~jls k aj )k (ou) are 
independent Af(0, s|sf)-laws, so (16.121) can be applied with (jj >k = s'j's k , 

E^^II^hsII^IIhs-^ 2 "- 

We also know that 

Si = maxcr^fc = ||S'||||S'|| ^ 1. 
From (I6.12p . we deduce that 

P{\\QA > a )< l Ch ~ 2n ~ ( 8 - 5 ) 
for some constant C > 0. Let 

M = dh~ n , (8.6) 

for some G\ ^> 1. Then (18. 5p gives 

PdlQllls > M 2 ) < Cexp (-r 2 VC), (8.7) 

for some new constant C > 0. 

We also want to control the trace class norm of Q u , so we will use the assumption 
that one of fh and fh is integrable. Assume for instance that m is integrable. Then 



m 



1/2 



is square integrable and we can factorize S = S1S2, with Sj G Op(m 1//2 ) being 
Hilbert-Schmidt operators. Let us write 

Now recall that the composition of two Hilbert-Schmidt operators is of trace class 
and the corresponding trace class norm does not exceed the product of the Hilbert- 
Schmidt norms of the two factors. Knowing that ||Si||hs = 0(h~ n ' 2 ) and applying 

flEZTD to s 2 J2j,k a jA uJ )^l s ^ we § et 

PdlQjtr > M 3 / 2 ) < Cexp (-h- 2n /C). (8.8) 
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In the following we will restrict the attention to Q^s with 

||QJhs<M, ||Q,J tr < M 3 / 2 , (8.9) 

and we have just seen that the probability that this is the case is bounded from 
below by 1 - Ce~ h ~ 2n / C . 

We wish to study the eigenvalue distribution of 

P s = P - 5Q^, (8.10) 

when 5 > is sufficiently small. (The minus sign is for notational convenience only.) 
Recall from Section [3], that for z e O, 

P(z) = {P-z)- l {P-z) (8.11) 

is a trace class perturbation of the identity. We now introduce 

P s (z) = (P- z)-\P - bQ^ -z) = P{z) - 8{P - z)- l Q w . (8.12) 

The Grushin problem will be used to find lower bounds for | det Ps(z)\. First we 
derive an upper bound: We have with Ps(z) = Ps, P = P(z): 

P5P5 (8.13) 

= P*P - 5(P*(P - z)- l Q u + Ql(P* - z)- x P - 5QKP* - z)-\P - z)- x QJ) 
= P*P + 5R, 

where 

\\R\\na < C(\\Qj ns + 8\\Qj\\Qj ns )<CM, (8.14) 
\\R\W < C(||Q w || tr + <J||Q w ||||Q w ||tr) <CM 3 / 2 , 

provided that <J||Q W || < 0(1), as will follow from (KTT5j) . 

In Section H] we studied P*P + «x(a~ 1 P*P) for h <C a -C 1. This operator is 
> a if l[o,i] < Xi as we ma y assume. Now assume that 

5M<^h. (8.15) 

Then 

Q, 

P*P + a X (a- 1 P*P) + 5R>-, 
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and 

lndetP;P 5 < \ndet(P£P s + a X (—)) 

a 

= lndet(P*P + ax(^—^) + 5R) 

a 

= \ndet(P*P + a X (— -)) + f tr ((P*P + a X {— -) + tR)- l R)dt. 

a Jo a 

The integral is C(l)|||P|| tr = 0(l)<JAff/a and combining this with (TC291 (assum- 
ing now (14.211) ). we get 

3 

lndetP^P < — ?— ( | \n\p\ 2 dxd£+ 0{l)a"]n±) + 0(1)^1. (8.16) 

Here we choose a = C7i , C 3> 1 and we can drop the last remainder term if we 
assume that 

5M§ < /i 1+ ^ n In 1 5 < ^1+^/2 ln 1 / 817 n 
rt ri 

For n > 2 this follows from (18.151) . but for n = 1 it might be a stronger assumption 
depending on the value of k. Then 

h*\detP s \<-±—(jj \n\p\dxdZ + 0(l)h K \n^). (8.18) 

Still with C a < 1 we define P + , P_, P = P, £ = £ as in Section [51 Here 2 
is fixed, P = P(z). With P s = P 5 {z), we put 

P<5 Pt- \ r 2/ D n\ (~\N , r2/nn\ w /-<Ar 



p 6 = : L (R n ) x C JV -> L (R n ) x C JV . (8.19) 

Now ||5(P - z)- 1 Q ul \\ < C5M, and 

5M h 
^<^<1 

under the assumption (I8.15p . so P5 has the inverse 



(8.20) 
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of norm < OiXj^fa). Writing 

= (P - z)- x Q u , (8.21) 
we have the Neumann series expansion 

Ss = (| (8.22) 



For < t < 5 we have 



|>„det Pi = ° 

0(i)4=IIQ u |lt, 



0(l)4=Af^ 

/a 



so 

lndetn = lndet(P)+0(l)-=Mi 
Applying f)5. 13j) . we get 

In |det7>a | = , \. ( [ [ \n\p\dxd£ + 0(l)a K \na) + 0(1)-?=M% . (8.23) 
(2iih) n J J V" 

Again, under the strengthened assumption ( 18.171) . we get with a = Ch, C 3> 1, 

ln|detP a | = J j \n\p\dxdi + 0{l)h*hx^). (8.24) 

The idea to get a lower bound for In | det Pg\ with high probability is now to use 
(12~TTD . (12321) which gives 

In | det P s \ = In | det P 5 | + In | det (8.25) 

and to get a lower bound for In | det E s _ + \. 
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9 Lower bounds on the determinant 



We keep the assumptions formulated in the beginning of Section [HJ, in particular 
(18.11) . We restrict the attention to the case when ( 18.91) holds with M given by (18.61) . 
and recall that so is the case with probability > 1 — Ce~ h " ' 1° . The restrictions 
(I8.15p . (18.171) on 5 will be further strengthened below. 

Using a formula of the type (18. 25ft we shall show that for every z e Q, the 
determinant of Pg(z) is very likely not to be too small. For that we study the 
probability distribution of the random matrix E s _ + , and show that we are close 
enough to the Gaussian case to be able to apply the results of Section [7| to the 
determinant. Recall that we work under the assumption ( 18. 91) . which is fulfilled 
with probability > 1 — Ce~ h n l c . We want to study the map 

oo 

Q^E s _ + = E\ + ^E _(5QEy~ l 5QEl 

i 

where Q = (P - z)~ x Q and \\Rj\\m < 1- Here, we used that \\E%.\\, < l/y/a. 

We can rewrite this further as 

Et + = E _ + + -(V^E°_Q^E° + + -= V -r fl j+2 9.1 





a 



We strengthen fl57T51> . (181?]) to 

SM 2 



« 1, (9.2) 



and recall that by ASM . M = dh- n . 
Then 

5Ch- n h n 

<-=■<!, 



^ C 

and we get 



— C 2 h~ 2n h 
Q = ^E°_Q^E° + + T, ||T||hs < ^— < 1. (9.3) 

'a 
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In view of (18.11) we have 

y^cE°_Q^El = ^E°_(P - z^S^W^JVaEl, (9.4) 

where we recall from (15.81) that 

N 

y^E + v + = J2Mj>j, sfaE°_v(j) = (v\fj), 1<]<N, (9.5) 
i 

where e%,..,epf and fi,...,f N are orthonormal bases for ran(l[ ^(P(z)*P(z))) and 
ran(l[o, a ](-P(^)-P(^)*)) respectively, writing ran(£?) for the range of B. 

Here, we wish to apply the discussion of Section [131 The operators SyfaE+, 
y/aE°_(P — z)~ 1 S are clearly Hilbert-Schmidt of rank < N. Let tj, tj denote the 
singular values of these operators so that tj — tj — for j > N + 1. 

Lemma 9.1 We have 

^<t v t 3 <C, 1 < j <N, (9.6) 
where C > is independent of h,a. 

Proof: (USD shows that 11^^+11, Hv^-ll < x > and clearly ||(P - ||5|| = 

Oil), so the upper bound in (19.61) is clear. 

On the other hand, ^/aE^ r v + is confined to a bounded region in phase space, 
and it is easy to show that 

C||SV^£>+|| > \\^E%v + \\ = \\v+\\, 

which implies that the smallest eigenvalue of ((Sy/aE^ )*(Sy/diE^_ )) 1 / 2 is > 1/C. 
The lower bound on tj follows. The argument for tj is essentially the same. □ 

Let /i, ...,/at and /i,...,/jv be orthonormal bases in of eigenfunctions of 
((^E°(P- zy^y/aE^P - z)- 1 ^)*) 1 / 2 and ((Sy^ElYiSy^El)) 1 / 2 respec- 
tively, with tj and tk as the corresponding eigenvalues. We can then choose the 
orthonormal bases {ej}, {e^} in L 2 so that 

e 3 = l(V^E°_(P - zy l Syf 3 , ej = l(S^E° + )fj, (9.7) 
for j = 1, 2, N. Then from (19.41) . we get 

V^E°_Q^E° + = 'j'^J-'-Ul- ( 9 - 8 ) 

l<j,k<N 
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Now we will be a a little more specific about the assumption (18. 9p . We will 
restrict the attention to the set Qm of matrices (atj,k{w)) such that 

\\S 2 Y, a j,k(u)e&S\\m < M, (9.9) 

which implies (I8.9P and which is fulfilled with probability > 1 — Cexp (— h~ 2n /C). 
Here we recall that we assumed m to be integrable and wrote S = SxS? with Sj G 
Op(S , (m 1/ ' 2 )). (When m is integrable instead, we make a corresponding factorization 
of S.) 

(19.31) can be reformulated as 

Q(a) = diag (tj) o no!j,fc)i<j,fe<jv + T(a.)j o diag (t k ), (9.10) 

||f(a.)||HS<0(l)^£, (9-11) 



. a 

for (a j>k ) E Q M - 

Let |||(aj,fc)||| denote the norm in (19.91) and let TC be the corresponding Hilbert 
space of N x N matrices. We shall view 138(0^, C^) =: TCn as a subspace of TC in 
the natural way. Note that the two norms are uniformly equivalent on this subspace. 

The Cauchy inequality implies (after decreasing M by a constant factor) that 
the differential of the map a. i— ► T(a.) satisfies the following estimate on Q M : 

\\df\\ n ^ nN = 0(1)^. (9.12) 



On TC, TCn we have the basic probability measures, 

m = ft (e-^ 2 ^^) m „ = n (e-l^l 2 ^^). (9.13) 

j,k=l j,k=l 

We shall now estimate II* (/z w ) on Q M , where 

n (("i,fe)) = ( a j,k)i<j,k<N + T(a.), (9.14) 

and to do so, we identify T(a.) with its image in TC under the natural inclusion 
TCn C TC, and write 

II = Il o k, n(a.) = a. + T(a.), U (a.) = (a j>k )i<j,k<N, (9.15) 
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for a. = (oij t k) G Ti. 

We first proceed formally, ignoring some technical difficulties due to the infinite 
dimension. We have 

III«HIIhs-H«-IIhsI = |||no«(a.)||m-||noa.|||a| (9.16) 

= | ||n K(a.)|| H s - ||n a.|| H s I x I ||n K(a.)|| H s + l|n a.|| H s I 

< ||f(a.)|| HS (2|||a.||| + 0(^)) 

< o(i) 6M3 



a 



'a 

where we strengthen the assumption ( 19.21) to 

SM 3 



< 1, or equivalents S < h 3n+1/2 . (9.17) 

As for the Jacobian of k, we recall that if A : Ti — > 7i is linear with ||A|| tr <C 1 
(uniformly with respect to M), then det(l + A) = 1 + 0(\\A\\ tT ). Also, if A is of 
rank < N 2 , we know that ||-A|| tr < iV 2 ||yl||, so in our case we get from (I9.12p that 

, dK JN 2 M 



Here the remainder term is <C 1 in view of the assumption (I9.17P and the fact that 
N < M. (Recall that N = 0{a R hr n ) by (142T1U 

If F is a locally defined holomorphic map : Ti — > Ti, then 

dF 

LidFix)) = I det — \ 2 L(dx), 
ox 

so in our case, 

5N 2 M 

Lidnix)) = (1 + _)L(dx). 

V« 

Combining this with (I9.16p . we get 

5M 3 

k*(hh) < (1 + 0(1)— p=-)^« on Q M . 



Since (JIo)*Hh — ^h n , we conclude that 



SM 3 

Mvn) < 1 + O l ^ ^, on Q M . 9.18 

'a 
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At the end of this section we shall complete the proof of (19 . 1 8[) by means of finite 
dimensional approximations. 

For a. = a.(u) G Qm we want to estimate the probability that |deti?i + | is 
small. According to Proposition I7.3[ the /^^(<iQ) -measure of the set of matrices Q 
with 

| det(diag (tj)- 1 o jE°_ + o diag (t,)" 1 + Q)\ < e a 

is 

< (7 e -|(C2+(iV-i)lnJV-2JV-a) 

if 

a<C 2 + (N + -)\nN-2N. (9.19) 

In view of (19. II) . ( 19. lOf) . ( 19.181) this is also (after a slight increase of C) an upper 
bound for the probability to have («,,&) G Qm and 

| det(diag (tj)- 1 o (jE°_ + + Q) o diag g)" 1 )! < e a , 



or equivalently that 

|det^ + |< e "(yn5II ? 



5 A A 
n 11' 1J' 

i i 



Write 

a = C 2 + (N - -)lnN -2N -b (9.20) 
and restrict the attention to b > 0. Then 

e a = e C2+(AT-i)lnAT-2JV-b 

and we get 

P( (P holds and |det££ + | < e JVini-ivini+(iv-|)i Il iv-Civ+c 2 -6) < e -6 (9>21) 

Summing up the discussion so far, we have 

Proposition 9.2 Consider the Grushin problem Ii8.19\) . Assume \^-21 ) and choose 
a = Ch, C ^> 0. Then there exist positive constants Co, C\, C*2, C such that for 
b>0 

Pf/TO) /ioWs and |det££ + | > e -5oh"- n hii-c 1 -Q a /»«-»i B 
> 1 -Ce- b -Ce- Coh - 2n . 



Here 5 is assumed to satisfy \9.1T\). 
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In view of f[Q4j) . <K25h . we get 

Theorem 9.3 We now return to the original P${z) in AS.lty) and we assume 
uniformly for all z in some open set Q <<= Q. If 5 satisfies Ii8.15\) . ( |<j?. there is a 
constant C > such that 

\n\detP s \<-^--(JJ \n\p\dxd£ + Ch K ln^), Vz G f2, (9.23) 

with probability > 1 — Ce~ c ° h 2 ™. If 5 satisfies the stronger condition (9.11), then 
there are constants C, C, Cq > such that for every z G Q and e > 0: 

In | det P s \ > — X —{ jj\n \p\dxd£ - Ch*(\n i + In ±) - e) (9.24) 

with probability > 1 — e - e ( 27rfc r™ - Ce~ c ° h ~ 2 " . 

Notice that the last term in the lower bound for the probability is much smaller 
than the second term, and can therefore be eliminated. 

We end this section by completing the proof of (19. 18)) by finite dimensional 
approximation. (We suggest the reader to proceed directly to Section [TUl) 

Lemma 9.4 We can choose the orthonormal bases {e^}, {e, } in L 2 so that \9. 7p is 
fulfilled for 1 < j < N and such that the square of the norm in h9.9\) is equivalent to 

J2\^M 2 Mj) 2 m 2 , (9-25) 

where faU), /^(/c) denote the singular values of S2 and S respectively. 

In this lemma we did not try to have any uniformity with respect to h. 
Assume the lemma for a while. Then for N > N + 1, we can replace in (18. ip 

by 

l<j,k<N jork>N+l 

which depends on finitely many random variables. Here a N (u) = (aj^uj))^ k< ^ 
and (3j k are the linear functions of a N which minimize 

\\S 2 o ( a hk^k + Yl PjJ&jtk) ^IIhs- 

j,k<N j or k>N 
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Here we can use the ej, ej of Lemma [9^41 On the set Qm, we have S 2 oQ N oS — > S 2 QS 
in Hilbert-Schmidt norm, and ||5 , 2Q jv 5'||hs < M when a(u) G Qm- 

We get the corresponding matrix = E°_ + + ^Q^ and can be written as 
in (I9.10p with T(a) replaced by T^(a) satisfying (19.111) . Now instead of /x^ we have 
the finite dimensional measure 



N 



Pn 



TT ( c -\a<. k \* L ( da j,k) ^ 
11 n 



N 



j,k=l 

which we can view as the restriction of /i^ to the tribe generated by ct^ with 
1 < j, k < N) and we define U N as in (19.141) with T replaced by T^. The subsequent 
arguments now become rigorous since we are in finite dimension and we get 

n?(A*K ff ) < (1 + 0(1)^=)ii Hn on Q M . 



Since II ^ — > U on Qm, we obtain (I9.18P in the limit. 

We next prove Lemma [9.41 
Proof: We consider first the following simplified problem. Let rris < 1 be a square 
integrable order function and let S G Op (S(ms)) be elliptic. We look for an or- 
thonormal basis e%, e 2 , ■ ■ in L 2 such that 



,u k Se k \\ 2 -^/i,(S) 2 K| 2 , (9.26) 
where fii(S) > p 2 (S) >...—> are the singular values of S and such that 



Since ftf = ®{h we nave Np 2 N = 0(h n ), and using also that fi N < 0(1), 



ei, eN 

is a prescribed orthonormal family of functions in S 

we get 

^ N - (Nh*y/ 2 + i' 

On the other hand there exists a constant kq > such that 

m s (p) > -^(p)- K0 , 

and we can use the mini-max principle to compare the eigenvalues of (S*S) 1 ^ 2 with 
those of (1 + ((hD) 2 + x 2 ))~ Ko/2 and deduce that 

1 1 
^ N ~ 0(1) {l + hN 1 '^/ 2 ' 
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If 0< ii <C 1, we have, with p denoting the symbol of (S*SY^ 2 , that 
dM(p, p -\[0,2n])) > dist^m^CtO,^])) 

> dist(0,{p; ^(p>-"° < |}) 



1 



If it is a corresponding normalized eigenfunction, we have (fi~ 1 (S*S) 1 ' 2 — 1)m = 0, 
and we notice that p:' 1 (S*S) 1 ' 2 G Op (S(p~ 1 ms)), where p~ 1 ms satisfies uniformly 
the axioms of an order function, when \i — > 0. We conclude that 

u = 0(1), in H(m) 

uniformly with respect to m if m — belongs to a family of orderfunctions that 
satisfy uniformly the axioms and m = 1 on {p G T*R"; p(p) < 2//}. From this, we 
deduce that 

fo>|u) = 0(ii N ), ViV, 

if G 5 is fixed, and [i — > 0. 

Let /i,/2, ... be an orthonormal basis of eigenfunctions of (S^S 1 ) 1 / 2 with /ii > 
/i 2 > ••• the corresponding decreasing enumeration of eigenvalues. Then (ej|/&) = 
0(k~°°), 1 < j < JV , k > 1. Let AT > JV . For j > JV + 1, put 

<& = - H^/, = /, + r„ E No 3 r, = (9(r°°). (9.27) 

Here, we let E No = (ei, ejv ) be the span of ei,...,ejv and 11^ : L 2 — > E No be 
the corresponding orthogonal projection. Then gj G and for j, k > N: 

(g J \g k )=5 hk + 0(r 00 k-°°). (9.28) 

Here the estimates are uniform with respect to A" and if A" is sufficiently large, we 
see that G = {{gj\gk)) is a positive definite matrix of which any real power has 
elements satisfying (BT2BJ) . Let (a j)k ) = G~ 1/2 , so that a j)k = 5 jik + O^' 00 ^ 00 ), 
j, k> N. Put 

e i = 3 > N - 

k>N 

Then e-,-, j > N form an orthonormal basis in the span of g^+i, gN+2, •••• We 
see that for j > N: 

e i =fi + J2 °U-°°k-°°)f k + r v E No 3 r 3 = 0{r°°)- (9-29) 

k>N 
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Gn = (Gjj) 1 - is a space of dimension N, containing E No . For 1 < j < N, we 
consider 

n Gjv /i = fj - n G i/j, 

oo 

JV+l 
oo 

JV+l 

oo oo 

JV+l £=JV+1 
oo 

= ^0(A;- 00 )/ fe + 0(iV— ), 

JV+l 

where the last term is in E No . Thus, we get for 1 < j < N: 

oo 

Kajj = fi + J2 °( k ~°°)h + ? v e n 3 r, = 0(N-°°). 

JV+l 

This implies 

oo 

KgJj =fj + J2 + N )~°°)f^ 1 < J < ^- (9-30) 
1 

Now, complete e\ , . . ., ejv to an orthonormal basis e\ , . . ., ejv in G^. Then e 1: e 2 , 

is an orthonormal basis in L 2 . (I9.30p shows that U.G N fi, ■■■,n GN f N is very close to 
being an orthonormal basis in Gn, and we see that 

JV oo 

e J = J2 u ^fk + J2°(( k + N y °)fk> 1<3<N, (9.31) 

k=l 1 

where (uj,k)i<j,k<N is a unitary matrix. (19.291) . (I9.3ip imply that for all j > 1: 

oo oo 

ei = E a i.*/* + E + jV )" DO ( A; + Ar )~°°)/^ ( 9 - 32 ) 

k=l k=l 

where a^k = u^k for j,k < N and a^k = 5j± when max(j, k) > N. 
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We now fix iV sufficiently large so that the above estimates hold. Using (19.321) . 
we get 

e; = fi + Y, °(r°°k-°°)h, Sej = njfj + ocr 00 * -00 )/*. 

k k 
(Se k \Sej) — i4;5j t k + O(k~ 00 j~ 00 ), 

00 oo oo oo 

1 k=l k=l 3=1 

= ((/i 2 + K)u\u)p 

= ((1 + K /i^ 1 ) fiu\fiu) 

where \i denotes the operator diag (//•,■). Here n~ l Kn~ l is compact: £ 2 — > £ 2 , so 
1 + ^T 1 K^,~ l is a non-negative selfadjoint Fredholm operator of index 0. If u(j) = 
0(j M °) and (l+n^K/i-^u = 0, then u = 0(J-°°). UO^vef, (1+ /i" 1 K n^v = 
0, we conclude that 

oo 

((1 + fi- 1 K f i- 1 )v\v) = ((/i 2 + K)u\u) = ||5(^ Mfe e fc )|| 2 > 

i 

thus 1 + fi^Kii" 1 is bijective and we finally conclude that (19. 26ft holds. 

Now we can finish the proof of the lemma. We choose {e?}, {ej} in (19.7p so that 
(I9.26P holds with S = S 2 , e,- = ej and 5 = S*, ej = ej respectively. Then the square 
of the norm (19.91) is equal to 

(S 2 ei\S 2 ej)(S*e k \S*e£,)ai ik aj t £ = (S <g> «Sa|a)^a, (9.33) 

where 

= (S 2 ei\S2ej), S £tk = (S*e k \S*e e ). 

From (I9.26P we know that 

S = /XPP/2, Jj, = diag Gu 2 (j)) 
S = J1PPJ1, Jl = diag Gu(j)), 

where P, P are positive selfadjoint operators satisfying 

— / <P,P < CI. 
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Then f !9.33|) can be written 

and the lemma follows. □ 

10 Spectral asymptotics when dp, dp are indepen- 
dent 

Let r <s Q be open with C 2 boundary and assume that for every z G dT: 

Tj z := p~ (z) is a smooth sub-manifold of T*R n on (10-1) 
which dp, dp are linearly independent at every point. 

This assumption, which is satisfied also in a neighborhood of dT, implies that 
codim (S 2 ) = 2. The assumption can also be rephrased more briefly by saying 
that dT does not contain any critical value of p : M. 2n — > M 2 . Here p is the leading 
symbol of the original (^-independent operator.) If p z {p) = (p(p) — z)^ 1 (p(p) — z) 
is the principal symbol of (P — z)~ 1 (P — z), we introduce 



I(z)= / \n\p z (p)\dp (10.2) 

which is the same integral as in (19. 23ft . (19.241) (where z was fixed). It is easy to 
see that I(z) is a smooth function on the neighborhood of dT where (llO.ip holds 
and as in [10] we can compute A z I(z). Since z i— > p z (p) is holomorphic, we know 
that A z ln |p 2 (p)| = when p z (p) ^ 0, ie when p £ S z . On the other hand p z (p) = 
{p{p) ~ — z ) where the first factor is holomorphic in z and non-vanishing, 

so 

A z In \p z {p)\ = A z In \p{p) — z\ =2ir5{z-p{p)). 
Uipe C^iVl), we get 

(A z I(z)Mz)L(dz) = J J A z (\n\ Pz (p)\Mz)L(dz)dp 
2n 5(z — p(p))(p(z)L(dz)dp = 2tt / ip(p(p))dp. 



Thus we get (as in [6l [7] when n — 1): 

1 



A{I{z))L{dz) =p*(dp) near dT, (10.3) 
2n 
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where dp is the symplectic volume element. Notice that this formula is still true 
without the assumption (110.11) and hence not only in a neighborhood of dT, but in 
Q; however I(z) is no more smooth in general but still well-defined as a distribution. 
This fact will be used in the proof of Theorem 110.11 

In view of ( 110. ip . we have V(t) ^ t and (I4.2ip holds uniformly with k — 1, when 
z varies in a neighborhood of dT. Correspondingly the conclusions in Theorem 19.31 
hold uniformly, when z varies in a small neighborhood of dT. 

Theorem 10.1 Let T d Q be open with C 2 boundary and make the assumption 
U0.1\) . Let 5 > satisfy ft9.17\) and assume that hhi^ <C e <C 1 (or equivalently 
5 > e~ e ^ Ch \ C > 1, e C 1, implying also that e > Chln^ for some C > 0). Then 
with C > large enough, the number N(Pg,T) of eigenvalues of Ps in T satisfies 

l^(n,r)-^voi(^(r))|<^ (io.4) 

with probability 

C _ 
> I g (2nh)« _ 

Proof: The eigenvalues of Ps in Q coincide with the zeros of the holomorphic 
function 

F s (z) =det P 5 (z). (10.5) 

Theorem 19.31 tells us that there exists a neighborhood of dil such that 
(a) With probability > 1 — Ce~ c ° h 2 ", we have 

In \F s (z)\ < + C/i In i), z G Q. 



(b) For every z e fl and e > we have 

hl|F ^ )l£ (2^ (/w - c '' (ln i +ln j ) - e) ' 

with probability > 1 - e - e(27rh)_ " - Ce~ c ° h ~ 2n . Notice here that In ± > In \. 

We can then repeat the arguments of (6j E]. Recall Proposition 6.1 from [7] 
proved in a more general form in [6]. 
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Proposition 10.2 Let Q, Q be open neighborhoods of dT and T respectively. Let 
ip G C°°(Q; R) and let f be a holomorphic function in Q such that 

\f{z;h)\ < e^ /% , zeQ,0<h<l. (10.6) 
Assume that for some e > 0, e C 1, 3z k 6 fi, 6 J, such that 

dTc\jD(z k ,V~e), #J = 0(-^), 

\f{z k -M>e> {Zk) - e \ keJ. (10.7) 
Then the number of zeros of f in T satisfies 

#(/- 1 (o) n r) = -L [ A V L(dz) + o{£), 

lnh Jr h 

where we let <p denote some distribution in V'(TU£l) extending the previous function 

The original statement in P, [7] was with a smooth function (p defined in a whole 
neighborhood of T satisfying (110. 6p there, but the proof works without any changes 
under the weaker assumptions above. 

In view of (a), (b), we can apply the proposition with h = (2irh) n and e replaced 
by 2e, ip = I(z) + Ch In ^, / = F$. Then ( 110. 6p holds with a probability as in (a), 
while ( 110. 7p holds with a probability 

>l-C-^h)-«_ Ce -c h-^_ 

We can define (p as a distribution in a full neighborhood of T by ( 110. 2p . Then 
' f AipL(dz) = — i— f ±-M{z)L(dz) = -\- 1 1 dxdi. 



2nhJv (2wh) n J r 2tt (2nh) n J J p - H r) 

The theorem follows. □ 

We next give a result about the simultaneous Weyl asymptotics for a family of 
domains 
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Theorem 10.3 Let Q be a family of domains T <<= f2 that satisfy the assumptions 
of Theorem \10.1\ uniformly in the following sense: Each T is of the form g(z) < 
(with g = gr) where g belongs to a bounded set in C 2 (Q) and g > 1/C on dQ and 
\dg\ > 1/C on dT, where C > is independent ofT. We also assume that U0.1\) 
holds for all z G dT, T G Q, uniformly with respect to (z, T). 
Choose 5,e as in Theorem \10.1\ Then with probability 

C e/a 

> 1 e ( 2irh ) n , 

e 



we have [TU^J) with a constant C independent ofT. 



Proof: As in the proof of Theorem llO.il we use Proposition ll0.3l with an appropriate 
grid of points Zk (see [7] for further details). We now need 0(1/ e) points to achieve 
that the union of the D(zk, \/e) covers the union of all the dT, rather than 0(1/ y/e) 
points as in the proof of Theorem 110.11 □ 



11 Counting zeros of holomorphic functions 

Let T (s C have smooth boundary dT. Assume for simplicity that 7 := dT is 
connected (or equivalently that T is simply connected). This is for notational con- 
venience only. For < r <C 1, we put 

7r = 1 + D(0,r) = dT + D(0,r). (11.1) 

Then 7 r has smooth boundary and is a thin domain of width w 2r. Let G r (z,w), 
P r (z,w) denote the Green and Poisson kernels of 7 r , so that the Dirichlet problem 

Au = v, u\ Qir = f, u,ve C°°(77), / G C°°(d lr ), 

has the unique solution 

u(z) — I G r (z,w)v(w)L(dw) + / P r (z,w)f(w)\dw\. 

We recall some properties of the Green kernel: If Q d C has a smooth boundary 
and Gq(x, y)is the corresponding Green kernel, then 

G n <0, (11.2) 
G n is C°° for x ^ y, (11.3) 
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Gq(x, y) = - - In \x — y\ for x y, X ->V & dQ } x — y — > 0. (H-4) 

x y 

£fi(~, -) = G rU (x,y), x,y e rft. (11.5) 

= -7,. is a very long domain of approximately constant width and (jll.4p is valid 
uniformly for x, y G Q, \x — y\ < 0(1), dist (x, 90), dist (y, dft) > 1/0(1). Moreover, 

\G r -^ r (x,y)\ < C e-^ Co , x,y E r' 1 ^, \x - y\ > (11.6) 

To recover these well-known facts, notice that r~ x ^ r is given by — 1 < (p(x) < 1, 
where (p(x) is the suitably signed distance from r~ 1 dT to x, so that |V<£>(:r)| = 1, 
|V 2 y3(x)| — 0(r). If u G ifo(r _1 7 r ), we have by integration by parts, 

((Vy?) 2 + ip(x)Aip)\u\ 2 dx = -2Re / y?(Vy2 • Vu)udx, 

implying 

1- 0(r))H 2 d£ < (2 + 0(r))||Vu|||H|, 
||u|| < (2 + 0(r))||Vu||, 
-A>(i-0(r)), 

where A = A r -i 7r is the Dirichlet Laplacian on r _1 7 r . From this estimate we can 
develop exponential decay estimates for —A, since we still have a positive lower 
bound for Re (e^(-A)e^) = - A - |W| 2 , if |V^(x)| 2 < 1/5. We drop the ensuing 
routine arguments. 

In view of (lll.5p . (111.61) we get 

\G r (x,y)\ <C e-^ lx - y [ x,ye 7r , \x - y\ > r/0(l), (11.7) 



\G r (x,y)\ = l±£&2ln|--H| > for \x-y\<r/C , (11.8) 
Z7r r r 

dist (x, 9Q), dist (y,dn) >r/C , |~ - ~| 0. 

Let (/? be a continuous subharmonic function defined in some neighborhood of 
TV- Let 

yU = n v = Ay? (11-9) 
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be the corresponding locally finite positive measure. 

Let u be a holomorphic function defined in a neighborhood of 77. We assume 
that 

h\n\u{z)\ < (p(z), z £ 7r- (11.10) 
Lemma 11.1 Let C\, C2 > 1 and let zq £ 7^ j_-> r 6e a point where 

hln\u(z )\ > (p(z ) - e, 0<e<l. (H-H) 
Then the number of zeros of u in D(z ,C 2 r) fl 7 (1 „j_y is 

< -f(e + J -G r (z , w)fi(dw)), (11.12) 
where C3 is independent of e, h. 

Proof: Writing (p as a uniform limit of an increasing sequence of smooth functions, 
we may assume that ip £ C°°. Let 

n u {dz) = S ^2n5{z — Zj), 

where Zj are the zeros of u counted with their multiplicity. We may assume that no 
Zj are situated on d^ r . Then, since A In \u\ = n u , 



h\n\u(z)\ = / G r (z,w)hn u (dw) + / P r (z,w)h\n\u(w)\\dw\ (11.13) 

' 7r J C*7r 



< / G r (z,w)hn u (dw) + / P r (z,w)tp(w)\dw\ 

J 7r •/ C?7r 



G r (z,w)hn u (dw) + (p(z) — / G r (z } w)n(dw). 

J It 

Putting z = z in (111.130 and using fill. Ill) , we get 

G r (z ,w)hn u (dw) < e + / —G r (z ,w)iJ l (dw). 

Now ^ 

-G r (z 0} w) > ——, C 3 > 0, 
27rC 3 

in D(^ , C 2 r) n 7(i_^) r > and we S et (PHI- D 
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Notice that this argument is basically the same as when using Jensen's formula to 
estimate the number of zeros of a holomorphic function in a disc. We could assume 
the bound Mn \u(z)\ < ip(z), in D(z , C 2 r) njnj_y =: Q r for some C 2 > C 2 . Then 

c 2 

we can replace the bound flll.l2p by 

C f 

T~( e+ L - G a.(^o,w)McM), 

which is sharper, since — Gq 1 < —Gn 2 , when Qi C il 2 . 
Now we sharpen the assumption (11 1.1 II) and assume 

hhy\u{zj)\ > (p(zj) - e, (11-14) 

where Zi, € 7(i„_j_) r are points such that 

N 1 
7, cljD^dr), N~-. (11.15) 

We may assume that z±, z 2 , z N are arranged in such a way that 

- z k \ ^ rdist (j, fc), j 7^ fc, (11.16) 

where j, are viewed as elements of Z/NZ and we take the natural distance on that 
set. We will also assume for a while that ip is smooth. 
According to Lemma 111.11 we have 

#( M - 1 (0) n {D( Zj , dr) n 7 ( i-^>)) < % e + / "OM*"))- ( 1L1? ) 

We introduce 

r=(l- — )r (11.18) 
w 

and consider the harmonic functions on 7f , 

#(*) = /i(ln|u(*)|+ / -Gr(«,w)nu(dw)), (11.19) 



7? 



$(2) = ¥>(*)+ / -Gr(z,w)ii(dw). (11.20) 



7? 



Then $(-2) > y?(-2) with equality on d^. Similarly, > /iln with equality 

on dj? . 
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Consider the harmonic function 

H(z) = $(z) - z G 7 f . (11.21) 

Then on dy?, we have by (jll.lOp that 

H(z) = <p(z) -h\n\u{z)\ > 0, 
so by the maximum principle, 

H(z) > 0, on 7f . (11.22) 

By f ill. 141) . we have 

H( Zj ) = <S>( Zj ) -¥( Zj ) (11.23) 
= <f(zj) — h\n \u(zj)\ + / —Gr{zj,w)n{dw)— I —Gr(zj,w)hn u (dw) 

J If A? 

< 6+ —Gr(Zj,w)lJ 1 (dw). 

Harnack's inequality implies that 

H{z) < 0(l)(e+ j -Gr{z j ,w)fi{dw)) on D{ Zj , dr) n 7(i-£)r- ( n - 24 ) 

Now assume that w extends to a holomorphic function in a neighborhood of 
T U77. We then would like to evaluate the number of zeros of u in V . Using (111.171) . 
we first have 

c N r 

#(^ 1 (0)n 7 f) < T (Ne + J2 / (-G r (z 3 ,w))fi(dw)). (11.25) 

i=i ^ 

Let x £ Co°(r U 7f 1 _i) ? r; [0, 1]) be equal to 1 on T. Of course x wu l have to 
depend on r but we may assume that for all fceN, and as r — > 0, 

V k X = 0{r~ k ). (11.26) 

We are interested in 

J x(z)hn u (dz) = J hln\u{z)\A x {z)L{dz), r = (1 - ^-)f. (11.27) 
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Here we have on 75? 

hln\u(z)\ = #(z)- -G 7 {z,w)hn u {dw) (11.28) 
= $(z) — H(z) — / -Gr(z,w)hn u (dw) 

= ip(z) + / —Gr(z,w)n(dw) — H(z) — I —Gr(z,w)hn u (dw) 

= <p(z)+R(z), 

where the last equality defines R(z). 
Inserting this in (jll.27p . we get 



x {z)hn u {dz) = J x{z)Kdz) + J R{z)A x {z)L{dz). (11.29) 

(Here we also used some extension of p to T with \x = A(p.) The task is now to 
estimate R(z) and the corresponding integral in (111.291) . Put 

Mj = Qj = D( Zj , Cir) n 7 r . (11.30) 

Using the exponential decay property (111.71) (equally valid for GV) we get for z G 
n 7 f , dist (z, d(D(zj, C x r) n > rfO(l): 

j -G ¥ (z,w)fi(dw) < I -Gr{z,w)fi{dw) + 0{l)^2M k e~^ lj ~ kl . (11.31) 

Similarly from (111.241) . we get 

H{z) <0(l)(e + / -G f (^,w)/u(^) + ^e~^ b " fc| M fc ), (11.32) 

for z G Oj n 7f . 

This gives the following estimate on the contribution from the first two terms in 
R{z) to the last integral in ([11.290 : 

I { I -Gr{z,w)n{dw) - H{z))A x {z)L{dz) (11.33) 
= 0(l)(Ne + Y^( sup I -Gr(z,w)fi(dw) + J2 e ~^ lj ~ klM k)) 



0(l)(iVe + ^ sup / -Gr(z,w)fi(dw) + /x(7 r )). 
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The contribution from the last term in R(z) (in (111.281) ) to the last integral in 
(TTT^j) is 

/ / G ¥ (z,w)hn u (dw)A X (z)L(dz). (11.34) 



Here 



G ¥ (z,w)(A X )(z)L(dz) 

Gr(rz, rw)A z x(r r z)r 2 L(dz) 



= j G?^(z,w)Az( x (rZ))L(dz) = 0(1), 

so the expression flll.34[) is 

O(/i)#(^ 1 (0)n 7 f) (11.35) 

N . 

= (!)(-+E / (-Griz^Mdw)) 
e N f 

= 0(1)(- + J2 -G r (zj,w)^(dw) + fi(^ r )). 

r 1=1 Jfy 



J 

Using all this in (lll.29p . we get 



x(z)hn u (dz) = J X (z)fJ,(dz) (11.36) 
+0(1)(- + Y^( sup / -Gr(z,w)fi(dw) + _-G r (zj,w)n(dw)) + fi(j r ))- 

We replace the smoothness assumption on (p by the assumption that tp is continuous 
near T and keep (11 1.14ft . Then by regularization, we still get (111.361) . 

In order to simplify this further, we introduce a weak regularity assumption 
on the measure /i. Assume first that fi = Aip is defined in a fixed r-independent 
neighborhood of dT. For D(z,t) contained in that neighborhood we assume that as 
t 0, 

W z (t):=fi(D(z,t)) = (D(tP°), (11.37) 

for some < po < 2. 
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Remark 11.2 It is easy to see that this assumption on Aip implies that ip is con- 
tinuous near T. In the case po > 1) we notice that as r — > 0, 

/i( 7r ) = O^ 1 ). (11.38) 

(This is true also for po < I but then of no interest.) 



Lemma 11.3 Assume ( 11.31 ) for some po g]0,2]. Then for every domain C 7 r 
and every z G Q fl 7( 1 _^) r; we have for < t < r/2: 

[ -G r {z, w)p(dw) < O{l)t p0 In - + 0(1) m(~)//(fi). (11.39) 
Jn t t 

Proof: Write 



—G r (z,w)p(dw) = / — G r (z, w)p(dw) + / —G r (z,w)p(dw). 

n JD(z,t)nn Jn\D(z,t) 

For \z — w\ > t, we have — G r (z, w) < 0(1) In ^ (cf ( I11.5P ). so the last integral is 
0(1) ln(0/i(fi). For w G D(z, t) n ft, we have 



-G r (z,w) < 0(1) In ■ ' 



2 — W 



hence 



D(z,t)nn 



-G r (z,w)p(dw) < 0(1) / ln-dW^(s) 



e»(i)(K-)w z ( s )]*+ f l -w z {s)ds) 

s Jo s 

0(l)t po ln^. 



□ 



Corollary 11.4 Under the same assumptions, we have for every N G N: 

/ -G r (z,w)fi(dw) < N (l)(r N + \n(-)p(n)). (11.40) 
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Proof: We just choose t = r M , < M G N and use that lnr~ M = Mlnr 1 □ 
If we assume (111.381) . then the corollary allows us to simplify (111.361) to 

J X {z)hn u {dz) = J X {zMdz) + 0(1)"- + N (r N + ln(i)/i( 7r )). (11.41) 

Summing up the discussion, we have proved 

Proposition 11.5 Let T <s C have smooth boundary and let (p be a continuous 
subharmonic function defined near T. Then we have the following result, valid uni- 
formly for < e < 1 ; < r < 1, < /i < 1: Let u be a holomorphic function, 
defined in T + D(0,r) with hln\u(z)\ < <p(z), z G dT + D(0,r) and assume that 
there exist Zi, Z;v G dT + D(0, |) such that 

N 1 
dT + D(0,r) C (J D(zj,2r), N - -, h\n \u(zj)\ > <p(zj) - e. (11.42) 

Then with Qj = D(zj, 2r) n (dT + D(0,r)), we have 

\#{u-\o) nr) - ^ f A<pL(dz)\ < 5ffl(l + M ( 7r ) (n.43) 

+ >( sup / — Gr.(z, w)fj l (dw) + / — G r (zj, w)iJ l (d'w) ). 

. «e%n(sr+£)(o,5)) Jn 3 -n(9r+D(o,§)) ~ in, ' 



//we assume also that (TTjFty holds for some < po < 2, then we have for every 
N > 0: 

|#(tr l (o) n r) - -L / a^l(^)| < (n.44) 

+ OjvCI)^ + in(-)//(ar + D(0, r)))). 
n r r 

Example 11.6 If <p is of class C 2 near the boundary, then (II 1.37ft is satisfied with 
Po = 2 and n(8T + D(0,r)) = 0(r). We choose iV = 1 so that the right hand side 
of (111.441) becomes 

— r^(- + rln-). 
a r r 

If we choose r = y'e, we get 

^(tt-^O) fir)-^| A^(z)L(dz)| < ^v^ln^. 
In this case we loose a factor lne -1 compared to Proposition 6.1 in [7]. 
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12 Spectral asymptotics in a more general case 



Let r <s Q be open with C°° boundary. For 2 in a neighborhood of dT and < 
s, t -C 1, we put 

V;(t)=Vol{peR 2n ; \p(p)-z\ 2 <t}, W z (s) = V z (s 2 ). (12.1) 

Recall that in any bounded domain in phase space, the symbols |p 2 (p)| 2 = Q z (p) 
and \p(p) — z\ 2 are uniformly of the same order of magnitude. If we replace \p(p) — z\ 2 
by q z (p) in (112.11) . we get a new function V K new (t) such that 

V z {^) < V^(t) < V z {Ct). (12.2) 

For the purposes of this paper, we can therefore identify the two functions and resort 
to the second definition whenever we find it convenient. Also, when z is fixed, we 
will sometimes write V(t) instead of V z (t). 
Our weak assumption, replacing (110.11) . is 

3k g]0, 1], such that V z {t) = 0(t K ), uniformly for z G neigh (dT), < t < 1. 

(12.3) 

Example 12.1 When (fTTHI) holds for z G neigh (dT), it is clear that (TT2"T3|) is 
fulfilled with k = 1 and in particular this is the case when {p, p} 7^ on p _1 (<9r). If 
z G <9£ \ then (110.11) cannot hold, so if dT fl 7^ 0, the best we can hope for 
is that 

\/z G p~\dT), either {p,p} ^ or {p, {p,p}} ^ 0. (12.4) 

This is the situation considered in the 1-dimensional case in [8] where determin- 
istic upper bounds on the density of the eigenvalues were obtained. Following some 
arguments there, we shall see that if fa!2.4\ ) holds, then U2.3\) holds with k — |. In 
fact, if we assume (112.41) and if p(po) = z G dT, we estimate the contribution to 
W Zo (t) from a neighborhood of po in the following way: 

If |{p,p}(po)| ^ t/C, then dRep, dlmp are independent near p and the contri- 
bution is 0(t 2 ). If |{p,p}(po)| is ver Y small, we know that \{p, {p, p}}(po)\ > 1/C 
and in order to fix the ideas we assume that H^ ep lmp(p ) > 1/C. This means that 

H:={p; {Rep,lmp}(p) = 0} 

is a smooth hypersurface in a neighborhood of po and that i?R e p is transversal to H 
there. A general point in a neighborhood of po can therefore we written 

p = exp (tH Rcp )(p), t G neigh (0, R), p G #. 
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Then Rep(p) = Kep(p'), lmp(p) = Imp(p') + t 2 g(t,p), g > 1/C. Write Rep = s so 
that a general point p' in H is parametrized by (s, p") with p" G neigh (0, R 2n ~ 2 ). 

Write z = xq + iy$. For every fixed p", if \p(p) — z Q \ < r then \s — xq\ < t and 
| Imp(p') + t 2 g(t, s, p") — 2/0 1 < r. Then we are confined to an interval of length 2r in 
the s-variable and, for every such fixed s, to an interval of length 0(t 1 / 2 ) in the t- 
variable, or to the union of two such intervals. By Fubini's theorem, the contribution 
to the volume is therefore C(r 3 / 2 ). Hence W zo (t) = C(r 3 / 2 ), so V zo (t) = C(t 3/4 ), as 
claimed. 

In Section [10] we introduced the distribution I(z) in (110.21) and showed that I(z) 
is subharmonic, satisfying (110. 3p . This implies that 

/ A(I(w))L(dw) = 27rW z (s), (12.5) 

JD(z,s) 

and (112.31) is equivalent to 

/ A(I(w))L(dw) = O(t po ) } uniformly for z G neigh (dT), < t < 1, (12.6) 

JD(z,t) 

with p = 2k g]0, 2]. This is precisely the condition (111.371) for / = <p, p = A J. In 
view of (112. 2p . the assumption (112.31) is also equivalent to requiring (14.211) to hold 
uniformly for z G neigh (dT) with q = q z — |(p — z)~ l (p — z)\ 2 . 
Consider the holomorphic function 

F 5 (z; h) = det P s (z), zett, (12.7) 

where we recall that Pg(z) = (P — z)^ x (P^ — z). Theorem 19.31 and its proof give: 

Proposition 12.2 Let 5 satisfy Then there exist constants C, C , C > 

such that 

(a) With probability > 1 — Ce~ c ° h 2 " , we have 

In \Fs(z; h)\ < j^(^) + Ch s ° In I), (12.8) 

for all z in some fixed neighborhood of dT. 

(b) For every z G neigh (dT), e > 0, we have 

In \F 5 (z; h)\ > j^(I(z) ~ Ch s °(\n I + In 1) - e), (12.9) 
with probability > 1 - Ce~ e[2 ^ n - Ce~ Coh ~ 2n . 
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For the upper bound (112. 8p . we recall that the upper bound (18.181) was obtained 
when ||Q||hs satisfies the estimate (18.91) and this event is independent of z. 

We can now apply Proposition lll.5[ with ip equal to / + Ch K In ~ and with h 
there replaced by (2nh) n , with e in fill. 14j) replaced by 

0(l)(h K \n^- + h K \n^ + e), 

and with e in (112.91) large enough, so that e is the dominant term in the last expres- 
sion. In other words, we take 

e>/i K ln^ (12.10) 
o 

using also that ln5 _1 > \nti~ 1 . 

For < r < 1, choose z±, and N as in the first part of (111.421) . Then in 
view of (b) in the proposition, the last estimate in (lll.42p (with h there replaced by 
(27rh) n ) holds for all j with a probability 

> i _ £ e -f (*rh)-» _ Ce- Coh ~ 2n . 
r 

The term 

1 r 

AipL(dz) 



2nh 

in ( 1 11.44ft becomes after the substitutions h i— > (2nh) n , ip i— > /: 

1 /" AI(z) _ _ 1 



(2vr/i) n 

where we also used (110.31) . 



Theorem 12.3 Let 5 satisfy [9J7\ ). Assume / TTOj) . with k g]0, 1]. Let N(P + 
5Qu, T) be the number of eigenvalues of P + dQ^ in V . Then for every fixed K > 
and for < r C 1: 

\N(P + Y) - -J— / / dx#| < (12.11) 
(27T/i) n JJ p -i(r) 

£(- + CV(r* + \n(-) [ [ dxdO), < r « 1, 

with probability 

> i _ £ e -f (2^)- (12 .i2) 
r 
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provided that 

/i K ln^<e<l, (12.13) 
o 

or equivalently, 

e -ch? <S, C> 1, e < 1, 
implying that e > Ch K \n j-, for some C > 0. 

In (112.111) we want the right hand side to be much smaller than h~ n so it is 
natural to assume that 

ln(-) [ [ dxdi = 0(r a "), r -> 0, (12.14) 

r J Jp- 1 (dr+D(0,r)) 

for some ao > 0. When k g]|,1], we automatically have (112. 14)) with any a £ 
]0, 2k — 1[. In the right hand side of (112.111) . we first choose N > ao and we choose 
r = e 1 /(i+«o) ) so t h a t e /r, r a ° = O^e 1 ^). Then the right hand side of (112. lip 
becomes 

C _^L- 

~~ h n 

So, if 1 ^> e > Ch K In i with C sufficiently large, and 5 is as in the theorem, then 

- ^ //_ 1(r) ^i s (12-15) 

with probability 

>\--^ e -h^)- n . (12.16) 

This expression is very close to 1 except possibly in the case k = 1, n = 1. In that 
case, we replace k by a strictly smaller value and choose 5, e as above. 

Theorem 12.4 Let Q be a family of domains T as in Theorem \12.3\ satifying the 
assumptions there uniformly ( cf Theorem \10.3\) and in particular we assume U2.3\) 
uniformly for all z in a neighborhood of the union of all the dT. Then we have 
M2.11\) with probability 

C 

> 1 - e -j2^hr 

provided that 



h K \n \ < e < 1. 
o 
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13 Appendix: Gaussian random variables in Hilbert 



spaces 



In this appendix we review some generalities about Gaussian random variables in 
Hilbert spaces that seem to be quite standard to probabilists. 

Let ofi, a 2 , ... be a sequence of independent J\f(0, l)-laws, and let 7i be a complex 
separable Hilbert space. 

Proposition 13.1 Let Vi,v 2 , ... <EH be a sequence of vectors such that \\ v j\\ 2 < 
oo, then if the sequence n x < n 2 < ... tends to oo sufficiently fast, we have that 



Let S(lj) denote the almost sure limit. If is another increasing sequence tending 
to infinity, such that the limit 



exists almost surely, then S(u) = S(u) a.s. 

Proof: Let (Q, P) be the underlying probability space. Then fj := aj(u)vj can be 
viewed as elements of L 2 (Q,7i) of norm \\vj\\. They are mutually orthogonal since 
the ctj are independent. We thus have an orthogonal sum a j v j which converges 
in L 2 (Q; Tt) and as usual, using the Chebyschev inequality, we deduce the existence 
of a sequence of partial sums that converges a.s. □ 

Let ei,e 2 ,.. and fi, f 2 , ■■■ be two orthonormal bases in 7i. Let ai(u),a 2 (uj), ... 
be independent complex A/"(0, l)-laws, and consider the formal vector cij(u))ej. 
Almost surely, {aj(u)}'^ is not in I 2 so our vector is not in 7i. However, 
then a.s., we can define the scalar product 




aij(u)vj exists almost surely (a.s.). 




i 





1 



as in Proposition 113.11 since {(ej\v)}f G £ 2 . 

We now look for random variables (3x{lo), (3 2 {u), ... such that 



oo 



oo 




(13.2) 



1 



1 



64 



in the sense that the formal scalar products with / 1( f 2 , ... are equal. This leads to 
the definition 



;i3.3) 



k=l 



which is well-defined as in Proposition I13.lt since k i— > {e^lf)) is in £ 2 . For every 
finite N, the variable 

JV 

5>*|/> fc (w) (13-4) 



fe=i 



has the density 



|c fe=i 



**l/;)| 2 



-M 2 /l(e fe |/i)| 2 



where * indicates convolution products. Hence the characteristic function (i.e. the 
Fourier transform) is 

1 N 

ex P (-7(EK c *l/i)l 2 )l^l 2 



fe=i 



so (I13.4p is a normal distribution Af(0, Y^k=i \( e k\fj)\ 2 )- The unitarity of the matrix 
((e&|/j)) then implies that is a A/"(0, l)-law. 

Proposition 13.2 (3j are independent N '(0,1) -laws. 



Proof: We have already seen that f3j are jV(0, l)-laws. To see that they are inde- 
pendent, we compute (using Proposition [13TT]) the joint distribution of fli,02, Pn- 
Write 

oo 



(iV) 



where 



( a hk \ 
\^N,kJ 



<Tj,k = ( e k\fj)- 
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ak(uj)i>k is a random variable with values in C and with the characteristic function 



7T 

1, 



-xp ( --INOI 2 



It follows that 



= exp | -- ^ ^ °z,k(?m,ki£ m ) • 

\ t=\ m=l J 

( 1 N N - \ 

\ 1=1 m=l / 

where <jj = (aj^'kLi £ But the <jj form an orthonormal system, so finally, 

XZ?<* k v k (0 =exp - i|^| 2 . 
This means that the joint distribution of /3i, ...,/3jy is 



and that /3i, ...,Pn are independent. □ 
The random variable (113. ip is an jV(0, ||t>|| 2 )-law. 

If f G 7i is any finite linear combination of the f), we know by construction that 



If v G W is arbitrary, we write v = v e + r e , where v e is a finite linear combination of 
the fj and ||r e || < e. We conclude that almost surely, 



J\ r e) 



Here the last two terms are jV(0, ||r e || )-laws and hence as small as we like with a 
probability as close as we like to 1, when e is small eneough. We conclude that 

oo oo 

(^aj{uj)ej\v) = (%2Pj{w)fj\v) a.s. (13.5) 
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Proposition 13.3 Let Ti, Ti be two separable Hilbert spaces and let T : Ti — > 7~t be 
a Hilbert- Schmidt operator. Let aj{uo)ej, (5j{u)fj be as above. Then T(^^° otj{uj)ej) 
is well defined a.s. and equal to T(J2T Pj( u )fj) a - s - 

Proof: We define T(J2T a A UJ ) e j) as a A UJ )'~^ e i m t ne sense 01 Proposition \T3A\ 
using that 

l^ll Te ill| = II t IIhs < °°. 
Notice also that for every v G TL we have a.s. 

00 oo 

( T (52 a 3( u ) e i)\ v ) = 5Z a jM( Te j» a - s - 

1 1 

oo 

= ^2<Xj(w){ej\T*v). 
i 

The same considerations apply to T(J2T ftj( u )fj) so * n y i ew °f Q13.5j) . for every 
v G H we have 

00 oo 

1 l 

We get the same conclusion a.s. simultaneously for all v in any countable set, and 
letting v = gi, g2, where gj form and orthonormal basis in Ti,, we conclude that 

00 oo 

1 1 

□ 

Now let £, J 7 , Q, Ti be separable Hilbert spaces and let T : £ — > J 7 , S : Q — > 7i be 
Hilbert-Schmidt operators. If / G .F, (7 g £/, we also denote by g, f the corresponding 
multiplication operators C 3 z 1-* zg, zf G so that /*it = Then gf*u = 

(u\f)g defines an operator : T — > G which has the Hilbert-Schmidt norm ||y||||/||. 
Let /j, j = 1, 2, gij, j = 1, 2, ... be orthonormal bases in J 7 , Q. Then {gjfk\Tk=x ls 
an orthonormal basis for the space HS(JF, Q)oi Hilbert-Schmidt operators T — >• 
Now, 

S gj f k T = (S gj )(T*f k y, 

and 

ii%/^ii 2 hs = n%irirAii 2 . 
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It follows that 

Eii^CTHs^ii^iiisimiis, 

and we conclude that 

HS(.F, Q) 3 A\-> SAT e HS(£, H) 
is a Hilbert-Schmidt operator. The earlier discussion can therefore be applied: 
Proposition 13.4 Let a^ k (uS) be independent Af (0,1) laws. Then 

j,k j,k 

is almost surely defined as a Hilbert-Schmidt operator. Moreover, ifgj, fk are new 
orthonormal bases in Q, T , then there exists a new set of independent M{0, l)-laws 
f3j,k(w) such that 

S o (^a 3 , fc ( W )^) o T = S o (EfaMgjfi) o T a.s. (13.6) 
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